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CHAPTER  I 
INTRODUCTION 

1.1.  Control  In  Manufacturing 

In  any  manu-f acturing  situation,  the  natural  variability 
o-f  a  process  results  in  parts  manufactured  with  variable 
dimensions.   This  variability  of  dimension  of    a  part  can  be 
expressed  in  terms  of  tolerance.   It  is  possible  to  reduce 
the  tolerance  to  a  minimum  by  applying  proper  control  over 
the  manufacturing  process.   But  the  cost  of  obtaining  a 
higher  degree  of  control  in  terms  of  smaller  tolerances 
mounts  rapidly  (see  Fig.  l.l)  C8D.   The  costs  can  also  be 
reduced   by  smart  design  with  manufacturing  in  mind,  the 
choice  of  an  especially  appropriate  machine,  an  efficient 
quality  control  system,  or  quantity  production  with  special 
machines  CBD. 

In  situations  where  complete  control  over  the  process  is 
not  possible,  or  extremely  prohibitive  in  costs,  it  is  a 
standard  practice   to  have  many  tolerance  groups  for  a  part 
depending  upon  the  clearance  requirements  of  the  mating  part. 
In  such  cases,  it  is  necessary  to  manufacture  or  order  the 
mating  components  with  corresponding  tolerances. 
Tolerance:  There  are  two  extreme  permissible  sizes  for  a 
dimension.   The  largest  possible  size  for  the  dimension  is 
known  as  the  upper  limit  whereas  the  smallest  permissible 


Tolerance 


Fig      1.1.     Graph    of    cost    of    production 
vs    tolerance    [8] 


size  is  called  the  lower  limit.   The  difference  between  the 
upper  and  lower  limit  of  a  dimension  is  called  the  tolerance. 

Tolerance  can  be  expressed  as  follows: 
Manimum  limit  =  Basic  size  +  upper  tolerance. 
Minimum  limit  =  Basic  size  +  lower  toleranci 


:e. 


1.2.  Alt»rn«tiv»  Component* 

Alternative  components  are  the  components  with  the  same 
basic  dimension  and  tolerance  but  with  different  upper  and 
lower  limits.   Main  component  is  one  of  the  alternative 
components  with  its  dimension  closest  to  the  basic  dimension 
and  whose  incidence  is  largest  among  all  the  alternative 
components. 

Exampl.:   Dimensions  of  a  40  mm.  diameter  shaft  made  on  a 
turning  machine  vary  from  38.7  to  40.3  mms.   This  shaft  mates 
with  a  bearing  whose  clearance  requirement  is  0.2  mms.   Now, 
there  are  3  tolerance  groups  of  the  same  shaft  namely,  grade 
l:  min.  dia  =  39.7  mms.,  max.  dia  =  39.9  mms.,  i.e..  40-0.3 
to  40-0.1  mms.,  grade  2:  40-0.1  to  40-^0.1  mms.,  grade  3: 
40.0.1  to  40.0.3  mms.  40  mm.  is  the  basic  dimension.   Grade 
1,   grade  2,  and  grade  3  form  the  alternative  components  for 

the  Shaft  with  40.0  mm  diameter.   The  main  component  has  the 

dimension  40-0.1  to  40+0.1  mms. 

Some  common  examples  of  alternative  components  are: 

pistons  with  the  same  basic  diameter  and  tolerance  but  with 

different  upper  and  lower  limits,  bearings  with  the  san 


ime 


inner  diameter  and  tolerance  but  with  dif-ferent  upper  and 
lower  limits. 

1.3.  Problvmm  With  AltarnAtlv*  componsntm 

Dealing  with  many  alternative  components  instead  of  one 
component  is  a  very  difficult  task,  especially  in  planning 
and  control  of  production.   Alternative  components  result  in 
wrong  planned  order  quantities  in  Materials  Requirement 
Planning  (MRP).   If  the  Bills  of  Materials  (BOM)  used  for 
calculating  Work-in-Process  (WIP)  quantities  is  the  same  as 
the  one  used  for  MRP,  then  WIP  figures  are    also  inaccurate. 
This  problem  is  discussed  in  more  detail  in  chapter  2. 

In  an  automobile  industry,  it  is  common  to  see  a  large 
number  of  components  each  having  several  alternative 
components.   But  such  a  large  number  makes  it  almost 
impossible  to  plan  and  control  the  production  of  these 
components.   In  cases  where  achieving  complete  control  over 
manufacturing  processes  is  prohibitive  in  costs  or 
impossible,  it  is  possible  to  minimize  the  production  of 
undesirable  alternative  components  and  thus  maximize  the 
production  of  main  components  by  the  proper  choice  of 
machines  available.   If  a  component  is  made  on  several 
similar  machines,  then  allocating  maximum  order  quantity  to 
machines  which  have  a  higher  degree  of  control  would  result 
in  a  larger  number  of  main  components  being  produced.   But 


given  a  machine,  it  is  extremely  di-f-ficult  to  determine  its 
degree  of  control  over  the  processes,  especially  in  the 
absence  o-f  past  data.   In  such  cases,  it  is  possible  to 
estimate  approximately  the  degree  of  control  (in  terms  o-f 
percentage  production  of  alternative  components)  based  on  the 
decision  maker's  experience  and  subjective  evaluation.   Under 
such  conditions,  it  is  also  di-f-ficult  to  quantify  goals  and 
constraints  precisely.   Such  an  imprecision  in  the  system 
gives  rise  to  a  situation  where  decisions  have  to  be  made  in 
a  fuzzy  environment.   Since  the  process  of  optimizing  the 
quantity  of  desirable  components  is  ongoing  and  the  goal  is 
to  be  achieved  over  the  length  of  a  number  of  time  periods, 
it  becomes  a  multi-stage  decision-making  process.   The  most 
common  tool  used  in  multi-stage  decision-making  is  Dynamic 
Programming  and  it  is  widely  used  in  our  study.    we  assume  a 
deterministic  system  under  control  to  start  with  and  solve 
the  problem  of  alternative  components  using  Dynamic 
Programming.   Later,  we  fuzzify  the  goals  and  constraints  to 
solve  the  same  problem.   In  the  next  step,  states  and 
controls  are  also  fuzzified  to  obtain  the  solution  for  a 
further  complicated  system. 

The  purpose  of  this  study  is  to  obtain  a  mathematical 
model  of  the  system  and  devise  suitable  methods  to  maximize 
the  production  quantity  of  the  main  component  assuming 

1.  a  deterministic  system  under  control 

2.  a  deterministic  system  under  control  in  a  fuzzy 


envi  ronment . 
3.   a  -fuzzy  system  under  control  in  a  -fuzzy  environment, 

1.4.  Overview  o-f  this  report 

The  next  section  in  this  chapter  deals  with  the 
literature  review  of  various  methods  available  for  multi- 
stage decision-making  with  varying  degrees  of  fuzziness. 
Application  of  these  methods  to  various  systems  are    also 
di  scussed. 

Chapter  2  discusses  the  nature  of  the  problem  in  detail. 
A  mathematical  model  of  the  problem  is  formulated  using 
dynamic  programming  in  a  deterministic  environment.   An 
example  of  the  problem  is  also  solved  using  the  model. 

Chapter  3  presents  a  basic  overview  of  fuzzy  set  theory. 
It  also  contains  all  the  basic  definitions,  operations,  and 
techniques  in  fuzzy  set  theory  which  are  used  in  this  study. 

In  chapter  4,  the  mechanism  of  multi-stage  decision- 
making is  presented.  An  example  of  the  problem  is  solved 
using  Bellman-Zadeh  model  assuming  the  crisp  system  under 
control  in  a  fuzzy  environment. 

In  chapter  5,  the  problem  is  considered  to  be  a  fuzzy 
system  under  control  in  a  fuzzy  environment.   Two  methods  of 
fuzzy  dynamic  programming  with  fuzzy  control  and  state 
variables  are  discussed  with  their  relative  merits  and  de- 
merits.  An  example  of  the  problem  is  also  solved  using  a 
modified  approach. 


Chapter  6  summarizes  the  study  and  the  results.   A 
discussion  on  the  results  follows. 

l.S.  Lit»ratur»  R»vi»w 

This  section  surveys  the  current  state  o-f  literature  on 
multi-stage  decision-making  in  a  fuzzy  environment.   Though 
the  problem  of  alternative  components  is  well  known  as  a 
planner's  nightmare  and  is  almost  unavoidable  in  a  precision 
manufacturing  system,  there  is  no  evidence  of  this  problem 
ever  being  studied  at  all.   There  is  no  literature  available 
on  the  effect  of  alternative  components  on  MRP  and  WIP  and 
their  inclusion  in  the  BOM. 

Since  its  introduction  in  1973  by  Joseph  Orlicky 
C21],  Materials  Requirement  Planning  (MRP)  has  replaced  all 
the  previously  existing  planning  methods.   Now,  it  is   a 
standard  and  most  effective  planning  tool  available  today. 
This  prompts  us  to  study  the  effect  of  alternative  components 
on  MRP. 

Since  the  publication  of  the  innovative  paper  "Fuzzy 
sets"  by  L.A.  Zadeh  in  1965,  more  than  4000  papers  have  been 
published  by  the  researchers  in  this  area.   Bellman  and  Zadeh 
C33  used  this  new  concept  of  fuzzy  sets  to  introduce  the 
process  of  decision-making  in  a  fuzzy  environment.   They 
contend  the  practice  of  solving  problems  involving 
imprecision  by  equating  imprecision  to  randomness.   To  quote 


a    ^ew    words    irotn    them    "Much    of    the    decision-making    in    the 
real    world    takes    place    in    an    environment    in    which    goals,     the 
constraints,     and    the    consequences    o-f    possible    actions    are    not 
known    precisely.       To    deal     quantitatively    with    imprecision,     we 
usually    employ    the    concepts    and    techniques    o-f    probability 
theory    and,     more    particularly,     the    tools    provided    by    the 
decision    theory,     control     theory    and    information    theory.        In 
doing    so,    we   are    tacitly    accepting    the   premise    that 
imprecision,    whatever    its    nature,    can    be   equated    with 
randomness.       This,     in    our    view,     is    a    questionable 
assumption" . 

"Specifically,    our    contention    is    that    there    is    a   need 
for    differentiation    between    randomness    and    fuzziness,    with 
the    latter    being    a    major    source    of    imprecision    in    many 
decision    processes...".       Bellman    and    Zadeh    formulated    a    fuzzy 
decision    model     in    which    a    fuzzy    decision    is    a    confluence    of 
fuzzy    goals    and    fuzzy    constraints.       They    also    developed    a 
multi-stage    decision-making    model     in    a    fuzzy    environment 
using    dynamic    programming.        This    model     is    used    as    a    basic 
tool     in    our    study.       Fuzzy    dynamic    programming    is    a    natural 
extension    of    dynamic    programming    in    conjunction    with    fuzzy 


sets. 


Bellman    and    Zadeh' s    model     provides    a    basic    framework    for 
multi-stage    decision-making    in    a    fuzzy    environment    where    the 
system    under    control     is    either    deterministic    or    stochastic 
with    implicitly    or    explicitly    defined    termination    time.       A 


B 


■further  study  on  stochastic  system  under  control  in  a  fuzzy 
environment  with  a  fuzzy  termination  time  was  done  by 
Kacprzyk  in  [8D. 

A  dynamic  programming  model  in  a  fuzzy  environment  where 
the  system  under  control  is  fuzzy,  is  presented  by  Baldwin 
and  Pilsworth  in  CID.   They  use  a  fuzzy  state  mapping  to 
denote  the  fuzziness  of  control  and  state  variables.   In 
order  to  overcome  the  high  dimensional  state,  they  use  a 
method  of  approximation  where  fuzzy  reference  sets  for  both 
control  state  are  used  by  making  use  of  conditional  fuzzy 

statements  of  the  IF THEN  form.   Kacprzyk  presents  a 

much  simpler  method  CIO]  using  branch-and-bound  algorithm 
with  state  and  control  as  fuzzy. 

Applications  of  fuzzy  dynamic  programming  are    novel  and 
diverse.   Esogbue  and  Bellman  C7:  studied  the  application  of 
fuzzy  dynamic  programming  to  a  wide  variety  of  problems  like 
the  budgeting  problem  in  cancer  research  [43,  R&D 
management  C5:,  health  care  delivery  problem  involving 
anesthesia  administration  problem  C7],  ill  defined  energy 
systems  C7:,  and  water  resource  system  and  cluster  analysis 
C7:.   Kacprzyk  and  Straszak  [113  apply  fuzzy  dynamic 
programming  in  formulating  the  problem  of  optimal  control  of 
regional  development  trajectories.   Kacprzyk  and  Staniewsky 
C12D  made  use  of  fuzzy  numbers  in  applying  dynamic 
programming  to  long-term  inventory  policy-making  in  a  fuzzy 


en 


vironment  to  obtain  optimal  inventory  policies.   Samoylenko 
C22:  used  the  fuEzy  dynamic  programming  for  optimization 
problems  associated  with  computer  network  design  and  routing. 
Moore  et  al .  C18D  used  the  fuzzy  dynamic  programming  in 
automatic  speech  recognition.   The  field  is  nevertheless 
still  embryonic. 
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CHAPTER  2 
ALTERNATIVE  COMPONENTS  IN  PLANNING  AND  CONTROL  -  A  MODEL 

2.1.  INTRODUCTION 

As  we  have  seen  in  section  1.1,  inability  to  achieve 
complete  control  over  manufacturing  processes  results  in 
alternative  components.   Alternative  components  pose  many 
problems  i n  an  industrial  environment.   Some  of  them  are: 
representation  in  the  Bills  of  Materials,  wrong  requirement 
quantities  in  MRP,  and  wrong  estimation  of  WIP.   The  next 
section  in  this  chapter  discusses  all  the  above  problems.   A 
problem  involving  alternative  components  encountered  in  an 
automobile  manufacturing  industry  is  presented.   A 
mathematical  model  of  this  problem  is  formulated  using 
dynamic  programming  approach.   An  example  of  this  problem  is 
also  solved  making  use  of  the  model. 

2.2  R»pr»s»ntation  of  •lt»rn«tiv»  components  in  BOM 

Since  alternative  components  are  also  inventory  items 
C21]  and  form  many  options  for  the  same  part,  it  is  essential 
to  represent  them  appropriately  in  BOM. 

Consider  the  example  given  in  section  1.1.   Suppose  this 
shaft  has  to  be  fitted  with  a  bearing  to  form  a  sub-assembly 
and  this  sub-assembly  has  to  be  used  in  an  assembly  with  a 
quantity-per  C213  of  1.   Then,  there  would  be  3  options  for 
the  bearing  since  bearing  has  3  alternative  shafts. 


11 


One  way  of  representing  alternative  components  in  BOM  is 
to  define  3  unique  assemblies  with  each  assembly  having  each 
of  these  three  alternative  shafts  for  the  bearing.   If  we 
follow  the  same  logic,  there  would  be  3  options  for  the  end- 
product.   If  there  is  a  large  number  of  parts  each  of  which 
having  several  alternative  components,  or  if  each  of  the 
alternative  sub-assemblies  have,  in  turn,  many  parts  having 
many  alternative  components,  then,  the  number  of  options  for 
the  end-product  would  be  enormously  high. 

Another  way  is  to  represent  the  main  component  with  the 
usual  quantity-per  and  the  rest  of  the  alternative  components 
with  quantity-per  as  zero.   This  will  ensure  that  all  the 
inventory  items  are  represented  in  BOM.   But  while 
calculating  the  requirements  via  MRP  logic,  there  would  be 
requirements  for  only  the  main  component  and  the  rest  of  the 
alternative  components  would  have  a  requirement  of  zero. 

The  standard  practice  is  to  represent  all  the 
alternative  components  under  the  same  assembly  with  the 
quantity-per  being  split  among  the  alternative  components  in 
proportion  with  their  past  usage  (see  Fig.  2.1  and  2.2). 
With  this  arrangement,  it  is  possible  to  avoid  numerous  end- 
products.   This  would  enable  to  calculate  the  requirements  of 
all  the  alternative  components  for  the  assembly.   This 
practice  is,  however,  not  fool-proof  and  results  in  wrong  MRP 
and  WIP  figures  when  the  degree  of  control  over  manufacturing 
processes  for  the  part  is  not  predictable. 


12 


Shaft    Assembly 


Bearinq    CI) 


O.  3 


Shiaft    Assiy.  1 
C grade    ID 


Shaf  t^Ar.sembl  y 


Shaft    Assiy.  2 
C grade    3) 


Shaft    CCID       Bearing      Shaft    CC2:>       Beari 


O.  P. 


Sl-iaf t    Assiy.  3 
C grade    3D 


CID 


CID 


CID 


ng   Shaft  CC3D   Bear! 


CID      CID 


ng 


CID 


Fig.  a.  a.  Representation  of  Alternati 


ve  Components  in  BOM 


13 


2.3  Wrong  R»qulr»m»nts  in  MRP 

Let  the  requirement  of  the  assembly  be  1000  units.   Let 
the  inventory  levels  of  CI,  C2,  and  C3  be  400,  100,  and  400 
numbers  respectively.   Requirements  of  alternative  components 
are    calculated  as  below; 

Gross  Requirement  =  Net  Requirement  of  the  sub-assembly 
*  Quantity-per  of  the  component 

Net  Requirement  =  Gross  Requirement  -  Stock-on-hand 
CI:   Gross  Requirement  =  1000  ♦  0.3  =  300 
Net  Requirement  =  300  -  400  =  -100 

=  0,  since  we  already 
have  100  in  excess. 

C2:   Gross  Requirement  =  1000  *  0.5  =  500 

Net  Requirement  =  500  -  100  =  400 

To  be  ordered  =  400 
C3:   Gross  Requirement  =  1000  *  0.2  =  200 

Net  Requirement  =  200  -  400  =  -200 

=  O  since  we  already 
have  200  in  excess. 

In  overall,  to  make  1000  assemblies,  we  needed  1000 

numbers  of  shafts.   We  already  have  900  on  hand.   So,  we 

needed  to  order  only  100  numbers.   But   we  ended  up  ordering 
400  ! 
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2.4.  Wrong  WIP  figurss 

I-f  the  consumption  o-f  CI,  C2,  and  C3  are  di-f-ferent  than 
the  assumed  quant i ty-per ,  then  WIP  quantities  of  alternative 
components  calculated  using  quantity-per  in  BOM  are 
i  naccurate. 

Let  the  number  o-f  sub-assemblies  made  be  1000.   Let  the 
actual  consumption  of  CI,  C2,  and  C3  be  400,  300,  and  300 
respectively.   I-f  we  have  already  issued  500  numbers  o-f  each 
component  to  the  assembly  shop,  then 
i.   Calculation  using  BOM 

CI:   WIP  =  500  -  1000  *  0.3  =  200 
C2:   WIP  =  500  -  1000  *  0.5  =  0 
C3:   WIP  =  500  -  1000  *  0.2  =  300 
ii.   Actual  computations 

CI:   WIP  =  500  -  400  =  100 
C2:   WIP  =  500  -  300  =  200 
C3:   WIP  =  500  -  300  =  200 
In  actuality,  WIP's  of  CI,  C2,  and  C3  are  100,  200,  and 
200  numbers  respectively.  But  according  to  the  calculations 
made  using  quantity-per  in  BOM,  these  figures  are  200,  0,  and 
200  respectively.   So,  the  WIP  figures  are  wrong! 

2.5.  Dynamic  Programming  Approach 

The  problems  due  to  alternative  components  discussed  in 
the  previous  section  are  due  to  the  variability  of  process  on 
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machines.   The   variability  could  be  due  to  some  intrinsic 
characteristics  of  the  machine,  or  due  to  the  skills  of  the 
operator  working  on  the  machines.   In  cases  where  the  parts 
are  made  on  many  similar  machines  for  the  purpose  of  mass 
production,  it  would  be  appropriate  to  assign  the  order 
quantities  to  the  machines  in  proportion  to  their  degrees  of 
control  so  that  the  overall  degree  of  control  is  manimized. 
This  would  obviously  result  in  the  maximum  production  of  the 
most  desirable  component,  thus  minimizing  the  effect  of 
alternative  components  on  MRP  and  WIP.   in  cases  where  the 
part  is  made  on  only  one  machine  in  many  shift  with  many 
operators  working  on  them,  assigning  the  order  quantities  to 
the  operators  in  proportion  with  their  degree  of  control  over 
the  machines  would  result  in  maximizing  the  quantity  of  the 
most  desirable  component.   The  order  quantities  assigned  to 
different  machines  or  different  operators  constitute  the 
alternative  plans.   Since  this  process  of  maximizing  is 
ongoing  and  takes  place  over  the  length  of  finite  number  of 
periods,  it  is  a  multi-stage  decision  process.   The 
alternative  plans  would  be  the  control  variables  and  the 
quantities  produced  of  alternative  components  during  each 
period  would  be  the  state  variables.   This  is  an  ideal 
situation  for  applying  dynamic  programming. 
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2.6.  Dynamic  Programming  Model  -  In  a  Crisp  Environment 

Consider  the  example  shown  in  section  1.1  and  the 
standard  BOM  representation  of  alternative  components  as  in 
section  2.  2.   Let  us  assume  that  the  shaf t.«^  are  manufactured 
on  2  similar  machines  for  the  purpose  of  mass  production,  each 
producing  its  own  mix  of  alternative  components. 

For  the  sake  of  simplicity,  let  us  assume  that  the 
capacity  of  both  the  machines  be  B  per  period,  where  I  is  the 
number  of  machines  =  1,2,  and  let  the  minimum  number  that  must 
be  produced  on  each  machine  be  b 
Assumpti  ons : 

1.  Total  demand  never  exceeds  total  capacity. 

2.  Demand  not  met  during  the  previous  period  is  carried  over 
to  the  next  period  unless  explicitly  stated. 

3.  Initial  inventory  is  zero. 

4.  Production  quantity  is  always  in  multiples  of  lOO. 
Let  Z  =  Objective  function  value  for  the  period  J. 

w^  =  weight  given  to  the  alternative  component  i  in  the 

objective  function. 
x  =  inventory  at  the  end  of  period  J. 
Qj,^^=  Total  quantity  produced  during  the  period  J  for  the 

alternative  k  on  the  machine  1. 
^tl  ~  production  percentage  of  alternative  component  i  on  the 

machine  1. 
q^^.,^  =  quantity  produced  during  period  j  of  the 
alternative   component  i  for  the  alternative  k 

'^^j    -   consumption  of  alternative  component  1  for  the  period  J. 
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D.  =  Total  demand  -for  the  part  -for  the  period  j. 

P.  =  Planned  requirement  for  the  period  j. 

0.  =  net  quantity  to  be  produced  during  the  period  j. 

I..,  =  Inventory  o-f  the  alternative  component  i  at  the  end  of 

the  period  j  -for  the  alternative  plan  k. 

Since  our  objective  is  to  maximize  the  quantity  of  the 
most  desirable  component  based  on  the  weights  o-f  desirability, 
we  have  the  objective  function 

n 

"^^  'i  \?^,   "v^iJ)c  <2.1) 

Now,  f (q.  .   )  =  max.  value  of  Z  from  stages 

j  ,  j  ■••1 ,  .  •  .  ,  n-1 ,  n    with    q.   .  as    the    quantity    produced    of 

v<j-i)k  '    ^ 

alternative  component  i  for  the  periods  beginning  from  j, 
j+1 , . . . ,  n-1 ,  n. 

Since  the  quantity  produced  of  each  alternative 
component  depends  upon  the  percentage  production  factor  for 
each  machine, 

%3c  =  V  V  <2.2) 

Inventory  remaining  for  the  alternative  component  i  at  the  end 
of  the  period  j  for  proposal  k  =  This  period's  production  + 
Previous  period's  inventory  -  this  period's  consumption  , 
i.e., 

^jk  =  %k  ■"  ^<j-i,k  -  ^ij       <2.3) 

Total  inventory  for  the  part  (inclusive  of  all  alternative 

n 

components)     =    k      =    V     I  (7    ai 

j         *'_      ijk  \^.~>  I 

^^%j-i>k^     "    ^^^    max...  max     [Z      +    Z         +...+    Z    3.        ...  (2.5) 
k  .         k  k  J  J+l  n 

J  J  +  l  n 
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Since  the  process  -flow  is  forward,  stage  i  is  calculated 
first  be-fore  calculating  stages  j +  1 ,  j +2,  .  .  .  ,  n.   Thus,  the 
results  of  stage  j  will  influence  the  results  of  stage  j+1, 
J+2,...,n,  but  not  vice-versa. 


Z.  +  max . . . max  iZ„  +...+  ZH 
j  +  1     n 

|z.  +. . .+  z  I  = 


.•.  f  ( q    .  )  =  max   Z  .  +  max  .  .  .  max  i  Z  .  _  + 

j  +  1     n 


but  max .  .  .  max  iZ .       +...+  Z   r  =  '^  .<P  >.) 

j  + 1     n 

.-.  f  (q     )  =  max  | Z .  +  f .  (q  ,  )   (2.6) 

J 

Now,  net  quantity  to  be  produced  during  the  period  j  = 

Quantity  planned  for  the  period  j  -  inventory  at  the  end  of 

the  period  j-1,  i.e., 

0.  =  P.  -  X  .  (2.7) 

J     J     J-* 

Q    =0,  -k(100)-b   (2.8) 

jki     J  1 

Q,   =0.-0,    (2.9) 

jk2      J      jkl 

The  equations  2.8  and  2.9  can  be  used  for  obtaining 

various  production  alternatives  for  each  period  which  is 

based  on  the  inventory  remaining  after  the  previous  period 

(x  )  and  the  planned  requirement  (P  )  for  the  current  period, 
j  J 


2.7.  Solution  Using  Dynamic  Programming  Model 

We  are    going  to  use  the  dynamic  programming  model 

derived  in  the  previous  section  to  solve  a  problem.   For  the 

sake  of  simplicity,  the  computations  will  be  made  over  a 

length  of  2  periods. 

Assume  a  production  plan  P  =  1000  units  for  the  first 

1 

20 


period  and  P^  =  900  +or  the  second  period.   Let  the  maximum 
production  capacities  o-f  the  machines  be  B  =  B   =  600  units. 
Let  the  minimum  batch  quantity  for  each  machine  b   =  b   =  400 
units.   Let  the  demand  for  period  1,  D  =  1000  units  with  the 
quantity  split  up  among  CI,  C2,  and  C3  as  ISX,  657.,  and  20*/. 
respectively.   D^  =  900  with  the  same  split  up.   Let  the 
percentage  production  factors  (s.^)  for  each  component  on  each 
machine  be  as  in  table  2.2. 

Now,  the  first  step  is  to  calculate  various  production 
alternatives.   This  is  done  by  using  2.8  and  2.9. 

P^  =  1000,  j  =  1. 
Plan  1:   Q^^^  =  lOOO  -  0(100)  -400  =  600 

°ii2  ^    ^^^^   ~    °jki  ^  1000-600  =  400. 
Plan  2:   Q^^^  =  lOOO  -  1(100)  -  400  =  500 
Q    =  500 

122 

Plan  3:   Q    =  400 

131 

Q    =  600 

132 

This  is  shown  in  table  2.3. 

Now,  let  us  assume  weights  for  the  alternative 
components.   Since  C2  is  the  most  desirable  one,  we  give  a 
weight  of  1.  Let  the  weights  for  CI  and  C3  be  -0.8  and  -0.6 
respectively.   The  negative  sign  indicates  the  undi si rabi 1 i ty . 
Then  the  objective  function  becomes 

max  Z  =  -O.Bq    +  q  ^  -  0.6q 
Table  2.4,  2.5,  and  2.6  show  the  detailed  calculations. 

Refering  to  table  2.6,  we  get 

max.  value  of  Z  =  807,   k*  =  {k°,  k *> ,  i.e. 
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Period     Qty.  to  be  assigned  to 
machine  1     machine  2 

1  400  300 

2  600  600 

Now,  suppose  we  define  some  constraints  and  goals  ior    the 
deterministic  system  as  follows: 
Goal:  Maximize  Z  >  750 
Constraints: 

1.  Shortage  of  each  alternative  components  must  be  less 
than  lOZ  of  demand. 

2.  Quantity  of  C2  manufactured  must  exceed  657.. 

3.  Overall  utilization  of  each  machine  must  be  >  85V.. 

The  alternative  (k  .k.)  satisfies  the  goal  and  the 

9   4  ^ 

constraints  2  &  3  but  the  shortage  of  C3  for  both  the  periods 
is  much  higher  than  lOX.   Looking  through  the  alternatives  in 
table  2.4  and  2.6,  we  do  not  find  any  alternatives  which 
satisfy  all  the  costraints  and  the  goal.   The  constraints 
seem  to  be   conflicting  each  other.   This  means  that  the 
decision  maker  does  not  have  any  choices  that  strictly  meet 
all  his  requirements.   So,  this  situation  does  not  provide  any 
rational  basis  to  make  decisions.   The  best  way  to  get  out  of 
this  deadlock  is  to  choose  the  alternative  which  satisfies 
the  goals  and  constraints  as  close  as  possible.   The  term  "as 
close  as  possible"  is  highly  subjective  and  indicates  that 
there  is  fuzziness  involved  in  the  environment  in  which  one 
has  to  make  a  decision.   Next  chapter  deals  with  decision 
making  in  a  fuzzy  environment. 
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Table   2.5.     Production    plan    alternatives    for    p*^riod   2 


Machi  nes 

Al ternati  vos 
1       2       3 

4 

1 
2 

600     500     400 
300     400     500 

300 

e.oo 

25 


!iii-     1^ 

1    >    -.   N 

N 

la  10        + 

K) 

1 

1                    0 

1^ 

III           N 

1  0)    L 

1  L    01 

1  a  a 

1  I.   ui        o_. 

10    3            N 

O 

IH-  0  •□ 

0- 

1        .-    0 

(0 

1 

IM   >  .-. 

1 

1  CP 

01 

1  c 

o 

O 

a< 

I-  c       7  to 

O   — 

t 

1  C    0    T3       U 

1 

CD    1 

■U 

\.~)  .» 

o 

U) 

1  <»   4J      1 

1  E  a     ^ 

o 

T    II 

L 

1  01    £        T 
1  L    3         i. 

o 

0 

in 

H. 

1      in      3  N 

1 

in  M 

>>   C     M         U 

o  a  1 

01 

1  t.    0 

N 

in 

•"• 

1  0    U     + 

in 

1    II 

n 

l-u 

t 

1  C    L         J< 

o 

H 

1  01   Hi        7 

N 

1  >  4J   cr 

+ 

in  in 

a> 

1    C    H-                    -■ 

o 

n  to 

c 

I"  m          u 

in 

1              ^ 

1    II 

E 

7 

E 

1           0- 

I-           >^ 

CM 

0 

c  Ul  » 

M 

L 

i"          II    ■* 

Q. 

1 

r       rj         1 

O 

u 

■•J                 1 

K) 

o 

— ^ 

n     1 

+ 

CD 

e 

?           i    U     1 

O 

<j 

0)           ^^              1 

in 

c 

u       cr           1 

II 

> 

3          .,                1 

Q 

■□       X           1 

0                           1 

o 

^          J         1 

•0 

a       lA            1 

(N 

o 

• 

N     1 

+ 

I^ 

tM 

>~           II        U       1 

O 

in 

'■-'                     1 

-0 

0) 

"■ax         1 

M 

II 

■u   0        ■'           1 

n 

C  ."   0-        ! 

t 

2  V.             1 

O 

o 

t- 

3    01              '^       1 

-0  in 

C3  a         (J     1 

+ 
o 

-< 

0- 

II 

L                               1 

0                          1 

*         -0                 1 

0                  1 

cm-            1 

D  .-    I.       ^          1 

-•  r  01      j<     1 

,^ 

Q.  -tj  a             1 

L  in                 1 

0    3           0          1 

*    0   T3         V        1 

-    0                  1 

C    >  .-                   1 

t    01    L                   1 

-^    L    01                 1 

a.  Q.  a.            1 

1        N 
I      «r 
1      1^ 

1        M       1 
1        -0       1 
1        hv       1 

1      O 
1      0- 
1      to 

1        O       1 
1        0-       1 

1      to     1 

1     m 

1        (N 

1      1 

1        O       1 
1        t       1 

1         1        1 

1     in 
1     •-• 
1      1 

1     in     1 

1      1     1 

1      o 
1     «• 

in     1 
<t     1 

1     f^ 
1     in 

1       M 

[N     1 
r^     1 
to     1 

1     in 

1        -0 

O       1 

in     1 

1      o 
1      o 
1     in 

O       1 
0-       1 

in     1 

1     in 
1     in 

O      1 
-0       1 

M       1 

«     1 

o 
in 


CD 

to 


in 
to 


o 
o 

-0 


in 

-0 


to 
tv 


-0 
0- 

to 


in 
I 


o 


o 

CD 


o 
If) 


1      to 
1     in 
1     r^ 

1        CD       1 
1       -0      1 
1      tv     1 

1        -0 
1        0- 
1      to 

1        -0       1 
1        0-       1 
1      to      1 

1      o 
1     •» 
1      1 

1     in    1 
1     in     1 

1     in 
1      1 

in    1 

1     in 
1      <s- 

O      1 

m     1 

1     r^ 
1     in 
1      to 

n     1 
r^     1 
to     1 

1     \n 

1        -0 

O      1 

in     1 

1      o 

1        (D 

1     in 

O       1 
0-       1 

in     1 

1     in     1 

1     in     1 

O      1 

fl    1 

1      ri     1 

to     1 

03 

to 


O 
O 
■0 


in 

0 


o 
t 


in 
I 


in 


M 

to 


O 


o 
in 


o 

UI 


1     t^ 
1      r^ 
1     rv 

1     r/     1 

1        0-       1 
t      r^      1 

1      o 

1        (M 

1     « 

1       O      1 
1      (1      1 
1      9      1 

1     in 
1     in 
1      1 

1       O      1 
1       1^      1 

I        1       1 

1     in 

t     in     1 

1      o 
1     in 

in     1 
in     1 

1     f-- 
1     in 
1      to 

n     1 

N       1 

to     1 

1     in 

1        -0 

O       1 

m     !• 

1      o 
1        0) 

1     in 

O       1 
U-       1 

in     1 

1     in 

1        Ul       1 

O       1 

n     1 

1      rj     1 

to      1 

o 
to 


o 


m 

IS 

I 


in 


o 

•0 


CD 
to 


o 
o 

-0 


in 
-n 


26 


CHAPTER  3 
FUNDAMENTALS  OF  FUZZY  SET  THEORY 

3.1.  Introduction 

So  -far  we  considered  the  situation  in  which  the  decision 

environment  and  the  parameters  were  known.   However  this  need 

not  be  always  the  case  and  some  sort  of  uncertainty  may  exist. 

Some  of  the  distinguished  cases  under  which  a  decision  is  made 

aret 

-Certainty 

-Risk 

-Uncertainty 

-Imprecision  (fuzziness) 
-Ignorance 
Risk  and  uncertainty  can  be  handled  by  the  application 

of  probability  theory.   In  case  of  imprecision,  the  state  can 

not  be  well  defined  because  of  lack  of  complete  knowledge  and 

the  vagueness  of  the  state.   For  example,  consider  the 

decision  situation  where  the  decision-maker  wants  to  keep  the 

inventory  levels  as  low  as  possible.   The  term  "  as  low  as 

possible"  is  a  vague  term  and  is  not  precise.   In  some  cases, 

it  is  impossible  to  assign  a  specific  value  to  this  term 

because  of  the  i mpreci seness  or  fuzziness  involved  in 

describing  the  state  inventory  levels.   Some  other  examples  of 

fuzzy  terms  are,     a  set  of  "tall"  persons,  a  set  of  "big"  cars 

etc.   In  such  cases  the  classical  probability  theory  can  not 

be  applied  because  it  can  handle  only  random  situations  rather 

than  vague  situations.   Moreover,  in  the  absence  of  past  data, 
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probability  theory  can  not  be  applied  as    there  would  be  no 
basis  on  which  probabilities  can  be  estimated.   Many  o-f  the 
situations  encountered  in  practice  come  under  this  category. 
The  major  breakthrough  came  when  L.A.  Zadeh  introduced  -fuzzy 
set  theory.   It  enabled  the  researchers  to  model 
mathematically  the  situations  involving  fuzziness.   Basically 
its  idea  is  as  follows:  The  conventional  set  may  be 
characterized  by  its  characteristic  function  which  takes  on 
values  only  0  and  1,  such  that  0  means  that  a  particular 
element  does  not  belong  to  the  set,  and  1  means  that  it 
belongs.   The  elements  may  there-fore  either  belong  or  not  to 
the  set,  i.e.  the  set  boundaries  are  crisp.   In  a  fuzzy  set, 
the  characteristic  function  is  replaced  by  the  membership 
function  which  takes  on  its  values  in  the  whole  interval 
CO,i:.   The  values  of  membership  function  express  the  degree 
to  which  elements   belong  to  the  fuzzy  set.   The  boundary  of  a 
fuzzy  set  is  therefore  not  clear  cut,  and  this  is  virtually 
what  we  need  to  characterize  fuzziness.   The  theory  of  fuzzy 
sets  provides  elegent,  adequate  and  simple  enough  tools  to 
accommodate  fuzziness  in  mathematical  modeling. 

This  chapter  discusses  the  fundamentals  of  fuzzy  set 
theory  and  the  basic  concepts  of  the  operators  and  tools  that 
are  used  throughout  this  study. 
3.2.  Some  Basic  Concepts  And  Definitions 
3.2. 1 .  Fuz  z  i  ness 

Fuzziness,  in  the  mathematical  sense,  is  a  description 
of  vagueness.   However,  fuzziness  should  not  be  mistaken  for 
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ambiguity,  uncertainty,  or  imprecision.   An  ambiguity  deals 
with  several  competing  descriptions  o-f  a  set  of  activities  or 
observations.   Uncertainty  is  the  use  o-f  random  variables  to 
represent  possible  values  o-f  activities  or  observations  in  a 
well  de-fined  set.   Imprecision  implies  the  existence  of  a 
sharp  tolerance  interval  for  the  description  of  a  set  of 
activities  or  observations.   Fuzziness  occurs  when  the 
physical  constraints  defining  the  tolerance  are  themselves 
elastic  or  ill-defined. 
3.2.2.  Fuzzy  Sets 

A  fuzzy  set  is  a  class  of  objects  (or  activities  or 
observations)  paired  with  their  respective  degrees  of 
membership  within  the  set. 

The  degrees  of  membership  usually  have  the  range  .C0,1] 
with  1  denoting  total  or  complete  membership. 

Mathematically,  if  X  is  a  collection  of  objects  denoted 
by  X,  then  a  fuzzy  set  A  =  {(x,  /j    (x))  xeX> 

(J    <x)  is  called  the  membership  function  or  degree  of 
membership  of  x  in  A. 

Fuzzy  sets  are  represented  in  two  ways.  One  way  is  to 
represent  them  in  the  form  of  A  =  (x,  u  (x)>.  The  other  way 
is  to  represent  them  as  A  =  iu    (x  ) /x  +    u    (x  ) /x  )  +...> 

A    1      1    '^A    2     2 

Exampl e: 

A  =  {<1,  0.1),   (2,  0.2),  (3,  0.6),  (4,  1)> 
A  =  <0.1/1  +  0.2/2  +  0.6/3  +  l/4> 
Both  representations  are  used  in  this  study. 
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3.2.3.  Basic  Properties  And  Operators 
1.  Properties 
a.  Containment 

A  -fuzzy  set  A  e  X  is  said  to  be  contained  in  a  ^uzzy  set 
B  c  X,  written  A  £  B,  i -f  and  only  if 

/j(x)<AJ(x),  Vx€X 

b.  Normality 

A    -fuzzy    set    A    c   X    is    said    to    be    normal     if    and    only    i -f 
max    ti    (x)    =    1 

X€X 

c.  Support 

The  support  of  a  fuzzy  set  A  c  X,  written  as  Sup  A  is 
defined  as    Sup  A  =  (x  e  X;  /l^  (x)  >  0> 


d.  Convexity 

A   fuzzy    set    A    is    convex    if 

jj    <Xx      +     (l-\>x    )     >   min    lu    <x    ),u    (x    )] 
A  ±  a        ^  '^A       1    '  "^A       a 

V  X  ,  X   e  X  and  \  e  C0,1]  where  min  denotes  minimum 


B.  a-Level  Set 

The  (crisp)  set  of  elements  which  belong  to  the  fuzzy 
set  A  at  least  to  the  degree  a   is  called  the  a-level  set. 

A      =    {x    €    X|  p    (x )     2   a> 

I 

A^   =    (x    e    X|/j    (X)     >    a>    is    called    "strong"    a-level    set, 

or    strong    a-cut. 
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Fig.    3.1.       A 


convex   furzy    set 


-I  X 


Fig.     3.2.       A   non-convex   fuzzy   set 
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2.  Operations 

As  in  the  conventional  set  theory,  the  basic  operations 
in  fuzzy  set  theory  are  the  complement,  union  and 
intersecti  on. 
a.  Complement 

The  complement  of  a  fuzzy  set  A  c^  X ,  written  -|A,  is 
defined  as  ij       (h)  =  1  -  /j  (h),  V  x  e  X 

-t  A  A 

see  Fig.  3.2. 

b.  Union 

The  union  of  two  fuzzy  sets,  A,B  c  X,  written  A  +  B,  is 
defined  as 

A-t-B  A       ^        B  ' 

where  "^"  is  the  maximum  operator 
Example: 

If  X  =  {1,2,3,4>,  and  A  =  0.2/1  +  0.5/2  +  O.B/3  +  1/4 
and 
B  =  1/1  +  0.8/2  +  0.5/3  +  0.2/4,  then 

A  +  B  =  1/1  +  0.8/2  +  0.8/3  +  1/4 

see  Fig.  3.3. 

c.  Intersection 

The  intersection  of  two  fuzzy  sets  A,  B  5  x,  written  A  n 
B,  is  defined  as 

AJ^^<x>  =  A^^<x)  ^  ^^^^),  V  X  e  X 

where  "^"    is  the  minimum  operator. 
Example: 

If  X  =  {1,2,3,4>  and  A  =  0.2/1  +  0.5/2  +  0.8/3  +  1/4  and 
B  =  1/1  +  0.8/2  +  0.5/3  +  0.2/4,  then 
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Fig.     3.3.     Complement 
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Fig.     3.4.     Uni 
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Fig.     3.5.     IriLoi  socLion 
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A  n  B  =  0.2/1  +  0.5/2  +  0.5/3  +  0.2/4 
see  Fig.  3.4. 

3.2.4.  Fuzzy  Numbers 

A  fuzzy  number  N  is  a  convex  normalized  set  o-f  the  real 
1  ine  [R  such  that 

1.  There  exists  exactly  one  x   e  (R  whose  membership 
function   li  (x  )  =  1. 

N    m 

2.  /J^<x^  is  piecewise  continuous 

where  x   is  the  mean  of  the  fuzzy  number. 
Example: 

The  following  fuzzy  set  is  a  fuzzy  number  A   =   (numbers 
close  to  5>  with  the  membership  function 


Ai^(x)  =   |, 


(-(x-5)^)/4  +  1,  3  <  X  <  5 


x/3  +  8/3  5  <  X  <  8 

0  el sewhere 


3.2.5.  Triangular  Fuzzy  Numbers 

The  triangular  fuzzy  number  is  derived  from  the  shape  of 
the  membership  function  graph.   The  graph  is  composed  of  2 
linear  functions  joined  at  the  mean  value  x  .   In  general,  a 
triangular  fuzzy  number  has  the  membership  function  of  the 
form 

(x-a  )/{a-a    )       a  <  x  <  a 

2  1  1  2 


i^(x)  =   I  (a^-x 


el sewhere 
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01       3345678 
Fig.     3.6.     Numbers    close?    to    3 


A 


2       10      12      3      4 


5      6 


Fig.     3.7.     Graph   of    the    trianqular     fuzzy 
number    < numbers   close   to   l> 
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Example: 

A  =  <  numbers  close  to  1>  with  a  membership  function  o-f 


^A< 


'•  0 


x/3  +  2/3, 
-x/4  +  5/4, 


-2   4x4.     1 
1  <  X  ^  5 

el sewhere 


3.2.6.  Fuzzy  Relations 

The  relation  is  one  o-f  the  key  concepts  in  mathematics. 
The  same  is  true  for  fuzzy  relation  in  fuzzy  mathematics. 

A  fuzzy  relation  R  between  two  non-fuzzy  sets  X  and  Y  is 
a  fuzzy  set  in  the  cartesian  product  X  x  Y,  i.e.  R  <::  X  x  Y. 
It  is  defined  as 


R  =  C  (x,y) ,  M.<x,y) / (x,y) > 


V  (x,y)  c  X  X  Y 


Example: 

Let  X  =  {Horse,  Donkey>  and  Y  =  {Mule,  Cow> . 
The  fuzzy  relation  R  labeled  "similarity"  may  be  as  follows 

R  =  "similarity"  =  0.8/ (Horse,  mule)  +  0.4/ (Horse,  cow)  + 
0.9/ (Donkey,  Mule)  +  0.2/ (Donkey,  Cow) 

The  particular  grades  of  membership  give  the  degree  to 
which  the  respective  animals  are  similar,  e.g.,  a  horse  and  a 
mule  are    similar  to  the  degree  O.B,  while  a  donkey  and  a  cow 
are    similar  to  the  degree  0.2  only.   It  can  be  represented  in 
the  matrix  form 


R  =  "similarity"  = 


horse 
donkey 


mul  e 

Cow 

0.8 

0.4 

0.9 

0.2 
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The    fuzzy    relation    is    evidently    a    -fuzzy    set    and    all     the 
■Fuzzy    set    properties,    operations,    etc.    are    applicable.       One   of 
the    most    important    property    is    the    max-min    composition. 

Consider    relations    R    c    X    x  Y    and    S    c   Y    x   Z. 
It    is    expressed    asR«Sc:XxZ,     such    that 

^Ros***'^'     ^    '"^^^     <M^<x»y)     '^mW.z))  Vx    eX,     z    eZ 

y€Y 

Example: 

If    X    =    {1,2),    Y    =    {1,2,3>,    and    Z    =    {1,2,3,4>,    then    the 
composition    of    two    fuzzy    relations    R    and    S    is 


R 

y 

1         2         3 
0.3    0.8      1    ' 

0.9    0.7    0.4 


1 

Y    2 

3 


S 

z 

12         3  4 

0.7    0.6    0.4  0.1 

0.4       1       0.7  0.2 

0.5    0.9    0.6  O.B 


R     .    S 

z 

13         3         4 
0.5    0.9    0.7    0.8 

0.7    0.7    0.7    0.4 


An   example   calculation    is    as   below. 
P^,g<2,3)    =         max  (u    (2,y)    ^^    <y,3)) 

y€{l,2,3>      "  * 

=     <P^<2,1)     Ap^(i,2))     ^    (Mj^<2,2)     A  ^^(2,3) 

<«    (2,3)     /N  M     <3,2)  ) 
=     (0.9    A   0.4)     ^    (0.7    A   0.7)     ^    (0.4    A   O.B) 
=    0.4    ^^   0.7    ^0.4 
=    0.7 
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3.2.7.  The    Extension    Principle 

The  extension  principle  is  one  o^  the  most  important  and 
powerful  tools  in  the  theory  o-f  fuzzy  set.   It  makes  it 
possible  to  extend  some  known  models  or  algorithms  involving 
non— fuzzy  variables  to  the  case  of    fuzzy  variables. 

'-^^  ^1  S  ^j»'--»'^   S  ^  ^^    some  fuzzy  sets  and  /: 
X  X.  .  .  xX  -»  Y  be  some  fuzzy  function.   Then,  due  to  the 
extension  principle,  the  fuzzy  set  B  c  Y  induced  by  fuzzy  sets 
A  ,...,A   through  /is 

n 

(J    <y)  =  max   A  u   (x  ) 
B  i=i  A   i 

I, 

(X  ,...,x  )  e  X  x,...,x   :  y  =  /(x  ....,x  ) 
Example: 

Let  x^=   {1,2,3>,   x^   =   {1,2,3,4},   /   represents   the 

addition,  i.e.  y  =  x  +  x     and 

1    z 

A^  =  <0.1/1  +  0.6/2  +  l/3>  and 
A^=  (0.6/1  +  1/2  +  0.5/3  +  0.1/4>,  then 
B  =  Aj+  A^  =  0.1/2  +  0.6/3  +  0.6/4  +  1/5  +  0.1/7 
This  principle  could  be  extended  to  other  operators  like   min 
and  max . 

3.2.8.  Fuzzy  Sets  of  Type  2 

A  fuzzy  set  of  type  2  is  defined  by  a  fuzzy  membership 
function,  the  grade  of  which  is  also  a  fuzzy  set  in  the  unit 
interval  CO, ID  rather  than  a  point  in  [0,11.   In  other  words, 
a  type  2  fuzzy  set  whose  membership  values  are    type  1  fuzzy 
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sets  on  CO, 1 D . 

The  fuzzy  membership  can  be  represented  as 

p  :  X  ■+  CO, ID   where  /j  (x)  is  the  -fuzzy  grade  being  a 

fuzzy  set  in  CO, ID  (or  in  the  subset  J  of  CO,  ID). 


3.2.9.  Operations  on  Type  2  Fuzzy  Sets 

The  extension  principle  can  be  used  to  define  the  basic 

operations  for  the  type  2  fuzzy  sets  C18]. 

Let 

A(x)  =  {(X,  M  (x))},  B(x)  =  ((x,  ij    (x))> 

where 

u    (x)  =  Hu  ,    u    (x))|x  €X;   u,M  (x)  e  CO, ID) 
"^A  V.'  '^u     '       '     v'  u. 

u  (x)  =  C(v,  u    <x))lx  eX:  v,  ^J    <x)  eCO.lD) 

The  u  and  v  are    degrees  of  membership  of  type  1  fuzzy  sets 
t       j 

and  u  (x)  and  u    (x)  are  their  membership  grades  respectively, 
'^u  V 

t  J 

1.  Union 

u  (x)  =u  <x)  Uij  (x) 

AUB  A  B 

=    Cw,  u  (w)  I  w    =    max     <u,     v>,    u,     v     eCO,lD> 

AUB  i.  J  V  j 

where 

u  (w)     =  Sup  min    iu    (x>,     u    (x)> 

AUB  u  V 

v  =  max  <  u.  ,v  >  x.  j 


2.     Intersection 

Lt        (x)    =   u    (x)    n  U    (x) 
aOb  ^a  ^a 

-    t(w,     (J    ^^<w)lw    =    min     (u,     v),     u,     v     e    C0,13> 

AlTB  i  j  i.  j 

where 
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''AriB^'^'   ^  ^"P         "^^^      t;L/^(H),  ^      (x)> 

w=mi  n<u  ,v  >  i       ^  ; 

«■   J  •' 

Type  2  -fuzzy  numbers  -follow  commutative  and  associative 
laws  and   follow  De  Morgan  laws.   However,  they  do  not  satis+y 
distributive,  absorption,  and  complementary  laws  C18D. 
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CHAPTER  4 
DYNAMIC  PROGRAMMING  APPROACH  -  IN  A  FUZZY  ENVIRONMENT 

4.1.  Introduction 

In  the  previous  chapter,  a    dynamic  programming  model  was 
formulated  and  a  solution  to  the  problem  of    alternative 
components  was  obtained  in  a  crisp  (deterministic) 
environment.   Most  often  the  real  life  situations  are  not 
deterministic.   Most  of  the  quantities  can  only  be  expressed 
in  terms  of  ranges  and  the  limits  of  these  ranges  are    not 
fixed  but  are  vaguely  expressed.   In  such  situations,  it  would 
be  most  appropriate  to  solve  the  problems  using  fuzzy  set 
theory.   This  chapter  deals  with  the  process  of 

decision-making  in  a  fuzzy  environment.   Bel Iman-Zadeh ' s  model 
of  dynamic  programming  in  a  fuzzy  environment  is  used  as  a 
basis  for  solving  the  problem.   A  branch-and-bound  algorithm 
proposed  by  Kacprzyk  is  also  dealt  with.   A  fuzzy  solution  is 
obtained  for  the  problem  of  alternative  components  using  the 
branch-and-bound  approach. 

4.2.  Decision-Making  in  a  Fuzzy  Environment 

A  general  conceptual  framework  for  decision-making  in  a 
fuzzy  environment  was  formulated  by  Bellman  and  Zadeh  C3]  in 
1970.   This  model  is  extensively  used  in  many  fields  due  to 
its  simplistic  approach. 
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4.2.1.  Fuzzy  Goal,  Fuzzy  Constraint,  and  Fuzzy  Decision 

A  fuzzy  environment  consists  o-f  -fuzzy  goals,  fuzzy 

constraints,  and  a  -fuzzy  decision. 

To  define  the  above  elements,  it  is  necessary  to 

introduce  some  set  of  possible  options  (or  choices, 

alternatives,  etc.)  denoted  by  X  =  {x>,  where  x  is  a  generic 

element  of  X.   The  elements  of  X  may  be  values  of 

expenditures,  alternative  courses  of  action,  etc.   X  contains 

all  possible,  or  all  relevent  options. 

Fuzzy  goal:   Fuzzy  goal  G  is  defined  as  a  convex  fuzzy  set  in 

the  set  of  options  X,  i.e.,  G  c  X  characterized  by  the 

membership  function  u    (x). 

a 

For  example,  a  fuzzy  goal  can  be  expressed  in  words  as 
"x  should  be  substantially  larger  than  10"  whose  membership 
function  is  subjectively  given  by 

p  (x)  =  0,    X  <  10, 

=  (1  +  (x  -  10)"^)"*,     X  >  10  (4.1) 

Fuzzy  constraint:  Fuzzy  constraint  C  is  defined  as  a  convex 

fuzzy  set  in  the  set  of  options  X,  i.e.,  C  c  X,  characterized 

by  membership  function  u    (x). 

c 

For  example,  a  fuzzy  constraint  could  be  "x  should  be  in 

the  vicinity  of  15"  which  can  be  represented  by  a  membership 

function 

^    (x)  =  (1  +  <x  -  5)*)"*  (4.2) 

c 

Now,  the  general  statement  of  the  decision-making 
problem  in  a  fuzzy  environment  is  "Attain  G  and  satisfy  C". 
Here,  "and"  corresponds  to  the  intersection  of  fuzzy  sets. 
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Fig.     4.1.     Fuzzy   goal:      "x   should    be 
substantially    larger     tlian    10". 


Fig.     4.2.     Fuzzy   Constraint:      "x    -.ho.ild    ho 
in    tlic^    vicinity    ol     1<J". 
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Fuzzy  decision:   Fuzzy  decision  can  be  de-fined  as  the  set 

resulting  from    the  intersection  of    G    and  C.   In  symbols, 

D  =  G  n  C  (4.3) 

or,  in  terms  o-f  membership  functions, 

IJ    ix)     =    ^    ix)    ^  ^    ix)  <4,4) 

D       a  c 

Note    that    -fuzzy    decision    D    is    convex    since    G    and    C    are 
convex.       The    operator     "^"    is    minimum    operator    and    the    decision 
is    min    type    -fuzzy    decision.       Min    type    fuzzy    decision 
represents    sa-f ety-f irst    approach. 

More  generally,  suppose  that  we  have  n  goals,  G  ,...,G  , 

1        n 

and  m  constraints,  C  C  ,  then  the  resultant  decision  is 

1         m 

the  intersection  of  the  n  goals  and  m  constraints.   That  is, 
D  =  G  n  G  n  ...  n  G  n  c^  n  c,  n  . . .  c      (4.5) 
and 

^^J^   =   t^a  ^   •  •  •    ^  ^^fj  ^  l^c   ^   '  ■  '    ^  tJ^  <4.6) 

1  n      1  m 

In  a  broader  sense, 
Decision  =  Con-fluence  o-f  goals  and  constraints 

In  the  con-fluence  diagram  in  Fig.  4.3,  notice  that  the 
highest  value  of  Pp<x)  could  be  less  than  1.   This  means  that 
there  is  no  option  which  fully  satisfies  both  the  fuzzy  goal 
and  fuzzy  constraint.   It  is  reasonable  in  many  instances  to 
choose  that  x  or  x's  which  have  maximal  grade  of  membership  in 
D.   Therefore,  the  optimal  (or  maximizing)  decision  is  defined 
as 

/J  (x  )  =  max  u    (x)  (4  7) 

xeX 
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4.2.2.  Deterministic  System  Under  Control 

Let  LIS  assume  that  the  system  under  control  is 
deterministic  and  its  temporal  evolution  is  given  by  the  state 
transition  equation 

H    =  /(x  ,  u  ),    t  =  O,  1,  ...  (4.8) 

where  x  ,  x    e  X  =  {x}  is  the  state  at  the  control  stage  t 
and  t+1  respectively,  and  u  e  U  =  {u>  is  the  control  at  time 
t.   The  initial  state  is  x   and  the  termination  time,  which  is 
fixed,  is  N  <  oo. 

4.2.3.  Multi-stage  Decision-Making  in  a  Fuzzy  Environment 

In  multi-stage  decision-making,  there  are  specified 

control  space  U  =  {u>  and  the  state  space  X  =  {x). 

In  the  beginning,  we  are  in  some  initial  control  u   e  U 

o 

which  is  subjected  to  a  fuzzy  constraint  u    (u  )  and  proceed  to 

CO 

a  fuzzy  state  x  e  X  through  some  known  cause  effect 

relationship  representing  the  system  under  control  S;  on  x  ,  a 

fuzzy  goal  /j  <u  )  and  attain  x   on  which  u  (x  )  is  imposed  and 
c   1  2  a      2 

so  on.   See  Fig.  4.4. 

Now,  to  obtain  a  set  of  optimum  alternatives 
(u  ,...,u   )   we  have 

O  N-l  ' 

A<„<u^,  •  .  .  ,u^,   )  =     max      (^^''(u  )  ^.  .  .  ^   fj^~  *^  lu         )     A 

D   O'      '   N-l  CO  c       N-l 

u  ,  .  .  .  ,  u 

O  N-l 

'^/j'*(/(x    ,u    )))  (4.9) 

a     N-l    N-l 

Notice  that  the  fuzzy  goal  G   is  imposed  on  the  final 
state.   The  last  two  right-hand-side  terms,  i.e.,  u*^~*(u   ) 

C       N-l 
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and    u    n(/<><         ,u        ))     depend    only    on    the    control     u  and    not 

a       ^         N-l'      N-l  ^  '  N-l 

on  the  other  (previous)  controls. 

The  maximization  in  4.9  may  be  split  into  two  parts: 
over  the  control  sequence  u  ,...,u    and  over  u   .i.e. 

O  N-2  N-l 

u(u,...,u   !><)=       max     CLi(u)^...^iJ    (u    )^ 

D    4  N-l'   O  CO  C       N-2 

u  ,  .  .  .  u 

O        N-2 

^  max  (u'*~*<u    )     ^n^ifix         ,u    )))D  (4.10) 

C       N-l        O       N-l     N-l 
U 
N-l 

N**  2 

Since  the  term  u  (u   )  depends  only  on  the  control 

C       N-2 

u   ,  then 

N-2' 

u    (u  ,...,u    Ix  )  =     max    Lu    (u    )     ^   .  .  .     ^  u  (u    )     ^ 

D    O  N-l  '  O  CO  C       N-3 

u  .....  u 

O  N-3 

^  max  (p'^~^<u   )  /N  max  (/j'^"*(u   )  ^iJ*(/U         ,  u   )))] 

C      N-2  C       N-l        a       N-l'    N-l 

u  u 

N-2  N-l 

(4.11) 
On  repeating  this  backward  iteration,  which  is  the 
characteristic  o-f  dynamic  programming,  the  -following 
recurrence  equations  are  obtained. 

N-  i  ^      V  N-  i  ,         >v    N-  i  -t-  1 

ia> 
u 


»-'-   (X   .)  =  max  {^  (u   .)  ^;j      (x     ))        (4.12) 

a       N-V  C       N-L        O  N-V.+1 


N-V 

and 


'^N-Ul  =  ^^^N-C   "n-1>  ^^-l^^ 


The  optimum  sequence  o-f  controls  to  be  determined  is 

given  by  successively  maximizing  values  o-f  u    in  (4.12). 

N-i 

The  above  model  o-f  dynamic  programming  recurrence 
equations  were  -formulated  by  Bellman  and  Zadeh  in  [3]  and  is 
often  known  as  Bel  1 man-Zadeh ' s  model  of  multi-stage 
decision-making  in  a  -fuzzy  environment. 
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4.2-4.  A  Simpler  Approach  For  Forward  Flow 
-  Branch-And-Bound  Method 

Like  any  other  dynamic  programming  models,  Bellman  and 
Zadeh's  model  suf-fers  from  the  "curse  o-f  di  menti  onal  i  ty" 
associated  with  dynamic  programming.   However,  the  model  is 
most  suitable  for  backward  flow  problems.   Here,  the  total 
number  oF  options  would  multiply  as  the  number  of  alternatives 
in  each  stage.   Therefore  it  is  not  suitable  for  solving 
problems  with  forward  flow,  i.e.  problems  that  progress  with 
time  as  it  would  be  cumbersome  to  handle  large  dimensions. 

Kacprzyk  CIO]  proposed  a  simpler  version  which  overcomes 

this  difficulty  and  is  known  as  branch-and-bound  method. 

The  procedure  starts  from  the  initial  state  x    .       A 

o 

control  u   is  applied  and  the  state  proceeds  to  «   and  so  on. 


1 


Finally,  the  state  ><   is  arrived  from  x    under  u      For  a 

N  N-1  N-1 

finite  control  space  this  process  may  be  represented  by  a 
decision  tree  (Fig.  4.5).   Let  us  consider  a  simple  example  to 
illustrate  this. 

Example:   If  U  =  ic^,     c^>  and  N  =  2,  then  the  decision  tree  is 
as  shown  in  Fig.  4.5.   The  nodes  are    associated  with  the 
particular  x.,  i  =  1,2  obtained  via  the  state  transition 
equation  4.8  while  the  edges  -with  values  of  controls  applied 
at  the  particular  control  stages. 

To  each  path  in  the  decision  tree,  i.e.  to  each  sequence 
of  edges  leading  from  x^  to  x^,  representing  the  consecutive 
controls  applied  and  states  attained,  there  corresponds  some 
value  of 

50 


Fig.     4.S.     Decision    troe    in    br  anch-.TncJ-h.ound 
method. 


1.  O 


Cll 


Fi-g.     4.6.     Reference?    graph 
for    shortage    of    CI 
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The  problem  is  now  to  find  an  optimum  sequence  of 
controls  u  ,u  ,  maximizing  (4.14).   Due  to  the  -finiteness  o-f 

the  control  space  assumed,  the  problem  reduces  to  -find 

id       Ik 
u  ,...,u     (in  the  general  case)  such  that 

O  N-l  ^ 

/J  (u   ,  .  .  .  ,  u   \x  )  =     max      (m°<u  )  ^ ^  i/''*^  (u   )  ^ 

D    O  N-l  '   O  CO  C       N-l 

u  .....  u 

O  N-l 

^   f/'lx     )  > 
a   N 

is  equivalent  to  finding  a  path  in  a  finite  decision  tree. 

The  solution  based  on  the  full  enumeration  of  all  possible 

paths  is  possible.   But  such  a  procedure  is  highly 

inefficient.   This  "blind"  traversal  of  the  decision  tree  is 

cumbersome  and  some  implicit  enumeration  technique  is  needed 

which  neglects  non-promising  paths  and  traverses  a  possibly  a 

small  set  of  promising  ones. 

Consider  the  decision  tree  as  shown  in  the  Fig.  4.5. 

Here,  we  traverse  the  decision  tree  from  x   through  the 

o       ^ 

consecutive  edges  (controls)  and  nodes  (states).   The  essence 

of  the  partial  enumeration  mechanism  is  basically  equivalent 

to  answering  the  following  question:  If  we  currently  arrive  at 

some  node  (state),  then  to  which  node  (out  of  those  traversed 

so  far)  should  we  most  rationally  add  next  edges  (controls)  to 

best  continue  the  tree  traversal?   Or  in  other  words,  what  is 

the  best  continuation  of  the  search  process  ? 

Let  us  first  denote 

V   =  /J  (u  ) 
o    '^c   o 
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^v  =  AJ^<"„)  ^•-     ■^U  (u)  =v  /"Ml^tu^)      (4.15) 

K        C     O  C     K  k-1      C     K 


V  =     Lt     iu     )     ^    ...     ^    u^~*(u    )   =  V     ^    U^~*(u    ) 
T-1     '^C    O  '^C       T-1        T-2     '^C       T-t 

V  =  p^Cu  )  -^  .  .  .  A  /j'^~*(u    )  A  ^'■(x  ) 
T     "^C    O  "^C      T-1      '^a        T 


=  V 
T 


-i^'^a^^^  =Ai„(u^,...,  u^_Jm^)         (4.16) 


Notice  that  if    we  consider  some  sequence  of  controls 

'^o'"''"k  ^  ^  ~  ^'  t*^^"  <^'-'®  to  an  evident  property  o-f  "A" 
operator,  we  have  -for  each  k  <  w  <  T  -  1 

^k  ^  ^„  =  ^k  '^  '^c***'^)c*i'  '^  •••  ^  ^c^"v>  (4.17) 

because  by  adding  to  v   any  -further  terms  through  "^"    we  can 

not  increase  the  value  of  v  .   We  also  have 

w 

Thus,  let  us  assume  that  we  are  at  k   control  stage  and  have 

traversed  so  far  some  nodes.   We  look  now  for  a  node,  among 

those  attained  so  far,  to  add  the  edges  (to  apply  controls). 

Obviously,  the  node  which  has  the  greatest  value  of  v  should 

i 

be  chosen  where  i  €  {l,...,k}.   The  other  ones  do  not  lead  to 
any  optimal  solution  because  for  any  next  controls  added,  they 
obviously  can  not  give  any  better  value  of  v  . 

2 

4.3.  Solution  for  Deterministic  System  Under  Control 

This  solution  approach  assumes  a  deterministic  system 
under  control  in  a  fuzzy  environment  with  fuzzy  goals  and 
fuzzy  constraints.   To  solve  the  problem  of  alternative 
components  already  described  earlier,  we  must  specify  the 
particular  fuzzy  constraints  and  the  objective  fuzzy  goals. 
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To  maintain  simplicity,  it  is  desirable  to  assume  the 
membership  -function  o-f  the  fuzzy  goals  and  constraints  in 
piecewise  linear  form.   In  this  form,  their  definition 
requires  the  specification  of  two  values  only:  the  aspiration 
level,  and  the  lowest  or  highest  possible  value.   Needless  to 
say,  these  values  are  the  result  of  subjective  evaluation  made 
by  the  experts  in  that  particular  area  or  the  decision-maker 
himself . 

Goal:   Maximize  Z  after  n  stages  (membership  grade  u^) 

a 

Constraints: 

1.  Shortage  (as  a  percentage  of  demand)  of  each 

alternative   component  must  be  as  low  as  possible  (u   ) 

civ 

2.  Quantity  of  the  main  component  C2  should  be  as  large 
as   possible  {u      ) . 

3.  Overall  machine  utilization  after  n  stages  for  all 
the   1  machines  should  be  as  high  as  possible  (u       ) 

C3l 

Reference  graphs  for  the  above  goals  and  constraints 
could  be  established  by  specifying  the  upper  and  lower  limits 
(aspiration  levels  and  highest/lowest  values).   In  the  graph, 
the  upper  limit  or  aspiration  level  has  the  highest  grade  of 
membehrship,  i.e.,  1,  and  the  lower  limit  has  the  grade  of 
membership  of  zero.   Grades  of  memberships  of  intermediate 
levels  of  goals  and  constraints  can  be  obtained  by 
interpolation  using  the  reference  graphs.   Assuming  a  linear 
relationship,  in  a  reference  graph,  the  upper  and  lower  limits 
are  joined  by  a  straight  line.   Fig.  4.6  to  Fig.  4.12  show  the 
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reference  graphs  -for  all  the  goals  and  constraints. 

Now,  ^    (u  ,...,u  )  =    max     (jj  ^  (j        ^  u       ^  u      )     ^ 

D    l'      '   N  Cll     C12     C13     C2 

U   .....  U 
1  N 

A  (p*  A  ^*   A  ^*   A  p*  )  A  .  .  .  A  (^^  A  pN   ^  ^   N     ^^_  ^^^ 
Cll     C12     C13     C2  C31     C3Z        O 

Notice  that  the  number  of  alternatives  in  each  stage 
depend  upon  the  production  plan  for  each  period.  Also  the 
constraint  of  overall  machine  utilization  and  the  goal  are 
imposed  after  N  stages. 

For  the  sake  of  simplicity,  we  shall  consider  only  two 

stages  as  the  dimension  of  the  problem  increases  with  the 

number  of  stages. 

Calcul ations: 

Stage  0,  u  =  1: 
o 

Constraint  1  (minimize  shortages  of  all  alternative 
components) : 

Looking  at  the  table  2.4,  we  get 
shortage  of  CI  =  0  =  0'/.  :*  u    =1  (from  reference  graph  4.6) 

Cll  ^    r- 

shortage  of  C2  =  10  =  1.5'/.  ■*   u°        =1  (from  refernce  graph  4.7) 

C12  ^   ^ 

shortage  of  C3  =  10  =  57.  -*  ^°       =1  (from  reference  graph  4.S) 

CIS  ^       ^ 

Constraint    2     (u      ) : 

C2 

Looking  at  the  table  2.4,  we  get 
quantity  of  C2  manufactured  =  640  =  647. 
Interpolating  using  the  reference  graph  4.9,  we  get 

Ai    =  0.6 

C2 

Similarly,  for  u   =  2,  we  have 

o 

^Cll  =   ^'   ^C12  =  ''      ^l.S     =  °-='   ^C2  =  °-"- 

for  u   =3,  we  have 


55 


1.  o 


C12 


Shortage  of  C2  -> 

Fig.  4.7.  Reference  graph  Jfor 
for  shortage  of  C2 


1.  O 


C13 


Shortage  of  C3  -> 

Fig.  4.8.  Reference  graph  for 
shortage  of  C3 


20% 
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1.0 


C2 


40% 
Qtjant.l  ty   of    C2    -* 


80;i 


Fig.     4.9.     Reforonce    graph    for     r^until  I  t  y 
of    C2    produced 


1.  0 


C31 


Utilization  of  machine  1  -♦ 

Fig.  4.10.  Reference  graph  for  utilization 
of  machiine  1 
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1.  o 


C32 


Utilization  of  machine  2 

Fig.  4.11.  Reference  graph  for  utiiizati 
of  machine  a 


on 
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1.  0 


M 


a 

O       COO  800 

Z  value  -» 
Fig.  4.12.  Reference  graph  for  Z  v,:»l..ie 


p°=0. 6 


Fig.     4.13.     Decision    tree    for     the    problem   of 
alternative    components 
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0  .         o  o  o 

Ai       =    1,     M  =    1,    p       =    O,    ij        =    0.65. 

Cll  Ci2  Ci9  C2 

Now,  using  the  equation  4. IB,  we  have 
p  (u  ,  u  ,  u  >  =    max    Cu  '^  u        ^   u        ^   u      1 

D        l'    2'    9  ^Cll    ^C12    '^C13    ^C2 

u  ,  u  ,  u 

1    2    3 

and  /j°(u  >  =  1  /N  1  A  1  /\  0.6  =  0.6 

D    1 

/J**(u  )  =  1  A  1  A  0.5  A  0.63  =  0.5 

D    2 

t/*  iu    )     =     1  A  1  A  0  A  0.65  =  0 

D    9 

max  tp  <u  ,u  ,u  ))  =  max  t0.6,0.5,0> 

=  0.6 
which  corresponds  to  u  . 

According  to  the  branch-and-bound  method,  the  most 
promising  alternative  is  u  .   So  we  traverse  along  this  path. 
Stage  =  1: 

Here,  we  assume  that  the  shortages  -from  the  current 
period  are  not  carried  forward  to  the  next  period.   We  assume 
that  the  orders  not  met  are    lost.   Re-fer  to  table  4.1. 
Assume  that  v  =    u        ^   u^      ^   u^        ^  u* 

C      Cll     C12     C13      C2 
v'(u  )  =  1  A  1  A  1  A  0.5B  =  0.58 

D    1 

vNu  )  =  1  A  1  A  0.78  ^   0.61  =  0.61 

D    2 

v*(u  >  =  1  A  1  A  0.22  ^  0.64 
D   a 

v'<u  )  =  1  A  1  A  o  A  0.67  =  O 

D    4 

Now,  /J  (u  ,u  ,u  ,u  )  =  max  tv*  '^  u*  ^  o*  ^  u^> 

D    1    2    3    4  D     "^Cai    '^C32    '^O 

Overall  machine  utilization  -factor  (for  u>   = 

1 

?l?[!_i°!I_ES!li°^_l_±_Ei^n_i°!I_E§!li°^_2 
Machine  capacity  for  2  periods 

Percentage  utilization  of  machine  1  =  -29_1_^29_  =  ^qOV. 

1  ,  o   1 

^(Z3±       t^*^!  ~  <from  reference  graphs  4.10  and  4.11) 

y.  utilization  for  machine  2  (for  u  )  =^92_1_522 —  =  507 

1       1200 
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u        (u.u)     =0.17. 

C32         1         1 


Simi lar ly, 


u         (u     , u    >  =     1  , 
C31         12 

u        (a    , u   )  =    1 , 

fj  (u     ,  u    )  =     1  , 

C34         14  ' 

also,     Z^    =  Z°    +    Z^ 


1  ,      O  1  ,  ^      ,^ 

u        (u    ,u   )     =    O. 42, 

C32         1  2 

u         (u    ,     u    )     =    0. 63 

C32         1  3 

u'      <u°,u*)     =    0.B3. 

C32         1         4 


Z    (u    ,u.)     =    732    .♦  u    (u    ,u   >     =    0.66     (  +  rom    the    reference    graph 
1       1  oil 

4. 12) . 

Simi 1 ar 1 y, 

£j*(u°,u')     =    0.74,       u*(Li°,u*)     =    O.Bl,       u*(u°,u*)     =    0.89. 
ai2  ai3  014 

/.   u*(u°,u*)    =    0.58    /^   1    /^  0.17    ^  0.66    =    0.66 

•^D        l'      1 

p*<u°,u*>    =    0.61    ^   1    ^  0.4    /N  0.74    =    0.42 

D         1         2 

M*<u'',u*)     =    0.22    ^    1    -^  0.63    /N  0.81    =    0.22 


/L(*(u°,u*)     =    0    A    1    A  o.e3    ^  0.89    =    0 

D         14 

max     {    0.17,    0.42,    0.22,    0>    =    0.42 

Cu     >     =     <U     ,      U    > 
1  2 
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CHAPTER  5 
FUZZY  SYSTEM  UNDER  CONTROL  IN  A  FUZZY  ENVIRONMENT 

5.1.  Introduction 

In  the  previous  chapter,  we  saw  that  in  the  case  of 
deterministic  system  under  control  the  initial  state,  final 
state  and  control  were  non-fuzzy  (crisp).   Therefore  we  could 
determine  the  grade  of  membership  in  the  fuzzy  goal  directly. 
However,  in  most  practical  situations,  both  the  state  and 
control  may  not  be  exactly  known.   In  such  cases,  the 
fuzziness  extends  to  the  control  and  state  variables  and  the 
resulting  system  is  called  the  fuzzy  system  under  control  in  a 
fuzzy  environment.   This  chapter  discusses  the  methods 
available  to  extend  Bel Iman-Zadeh ' s  model  to  such  systems. 
Two  methods  proposed  by  Baldwin  and  Pilsworth  [1]  and  Kacprzyk 
CIO]  are  described.   The  problem  of  alternative  components  is 
solved  using  the  latter  method  with  some  modifications. 

5.2.  Fuzzy  System  Under  Control 

In  fuzzy  system  under  control,  both  the  initial  state 
and  controls  are  fuzzy,  and  hence  the  final  state  obtained 
becomes  fuzzy.   So,  the  grade  of  membership  in  the  fuzzy  goal 
can  not  be  determined  in  a  direct  way.   Some  kind  of 
"trickery"  is  needed  to  circumvent  this  problem  in  order  to 
obtain  the  decision  in  terms  if  a  point  value  of  grade  of 
membership.   Both  the  methods  presented  here  involve  some  kind 
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of  trickery, 

5.2.1.  Method  of  Baldwin  and  Pilsworth 

This  method  extends  the  Bellman-Zadeh  model  to  include 
fuzzy  state  mapping  which  allows  the  system  dynamics  to  be 
described  imprecisely. 
Some  basic  definitions: 

1.  A  fuzzy  mapping  /:  X  ^  Y  is  a  fuzzy  set  on  X  x  Y  with 
membership  function  /j  <x,y). 

2.  A  fuzzy  function  /(k)  is  a  fuzzy  set  on  Y  with  membership 
function  AJj,^^j<y)  =  A^^<x,y). 

3.  The  fuzzy  set  /(A)  on  Y  i s  a  fuzzy  mapping  of  a  fuzzy  set 

defined  as  p^^^^(y)  =  ^ax  (min  <u^x),    u  (x,y)),  V  y  e  Y>. 

xeX         A       f 

Consider  the  finite  fuzzy  automation.   Given: 
<i).  J    =    <U,  X,  /,  F(U>,  F(X)>  with  final  state  x(T>,  where  U 
and  X  are  finite  spaces  for  the  decision  and  state  sets, 
respectively,  and  /:  X  x  U  ->  X  and  F(U),  F(X)  are  the  sets  for 
the  fuzzy  decision  and  states,  respectively, 
(ii)  a  fuzzy  dynamical  transition  equation: 
^x<t-M>^'*^^*^'>  =  ^f,k>  =  max  (min  ifuik)  , 

A/j<k,x(t  +  l)))>,  (5,1) 

V  x(t^l)  €  X  and  t  =  0,1,..., T-1  where  J  ^  X  x  U;  k  ^  K 
and  J  is  a  fuzzy  set  on  K  representing  fuzzy  state  x(t)  and 
fazzy  decision  u(t)  with  membership  function  m.(x,u)  = 
min  <A^^,^,(x).  P^^^,(u)). 

Also  consider  fuzzy  constraints  i^^,    p^,  .  .  .  ,  ^^_^>,  ^^  ^  ^  ,^^  ^ 
where  p.  is  imposed  on  the  input  u.,  0  S  i  ^  T-1  and  also  that 
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a  fuzzy  goal  fj       is  imposed  on  the  -final  state  x(T). 

Define  a  relational  matrix  R(u,x)  =  ^    (u)  ^  /j  (x); 

o        a 

u  €  U;  X  e  X.   This  can  be  thought  of  as  representing  how  well 
the  control  and  goal  constraints  are  satisfied.   Thus,  for  a 
fuzzy  control  u   and  resulting  state  x  ,  we  can  form 
T(u  R  X  >  =  u  o  R  o  X  (5.2) 

O     1         O  i 

where    <•    represents    max-min    composition    as    given    in    3.2.6.       The 

decision    problem    is    then    to    choose   u    e   U       (the    allowed    set    of 

fuzzy    decisions)    such    that    T(u   R    x   )     is    maximized.       Thus, 

o  1 

T(u    R    X    )     =    max     [{max     (^^(x)     ^  fu    (,u)     ^      n       (u))>]    ^  ^      (x) 


xeX         u^J  o  1 


(5.3) 


=    (l/<^^)^(,^o^)  (5.4) 

a  X  o         u 

1  o 

according    to    lemma    1    in    CID. 

Define    a    functional    f    (x    >    =    max     .    T(u u       R    x    ),  (5.5) 

Jc   k  .  ,*       k        T-1     T 

given  that  the  system  starts  at  stage  k  in  fuzzy  state  x   with 
X   the  fuzzy  terminal  state  resulting  from  a  sequential 

T 

application  of  fuzzy  decisions  u  ,...,u    under  a  fuzzy 

rr-  k  T-1 

mapping  /.   By  invocation  of  principle  of  optimality,  a 

general  recurrence  equation  is  obtained  as  below: 

for    k    =    T-1,     T-2,...,     1,     0: 

if    (X    )    =    max     ^Z(u   o    tj      )    ^  r        (x        )];  (5.6) 

^k      k  _,,*         k         u  ^K+i      k+i       ' 

"k^k  " 

and    for    k    =    T,     the    terminal     stage: 

if     (x    )     =     (p^    *    ^      )  (5.7) 

T      T  a  X 

T 

Note  that  for  fuzzy  subset  x    of  (5.6)  the  membership 

k*l 
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."^ 


function    u  («        )     is    given    by: 

k+i 

l:        (x       >    =         max  (;::      (X    )    '^  ;.'      (u.  )    ^    (^'(x,,u,;x,      )) 

•^X  k+l  V    II    X,        It  Ll,        ^  f       k'     k'     k+1 

k+1  x.ueXxUk  k 

k        k 

=    fj      (x,)    ^  M      (u.  )     •     (x,<x,,u,;x,      ))  <5.e) 

X.        k  u,        k  f       k'     k'     k+i 

k  K 

It  can  be  seen  that  dynamic  programming  functional  equation 

for  this  problem  with  fuzzy  mapping  and  initial  state  in 

addition  to  fuzzy  environment  has  associated  with  it  the  well 

known  "curse  of  dimensionality".   To  overcome  this  difficulty, 

Baldwin  and  Pilsworth  propose  a  process  called  "fuzzy 

interpolation"  CID.   This  is  to  find  a  relational  matrix  for 

f   (x   )  using  such  statements  as  "if  x.    is  large  then 
k+i   k+i  k+i 

?   (x   )  is  small",  etc.  where  large  and  small  are  fuzzy 
k+i   k+i 

linguistic  variables.   With  this,  a  fuzzy  value  of  f   (x   ) 

k+l   k*l 

can  be  determined  for  any  given  x   .   However,  the  method 

k+l 

becomes  very  cumbersome  and  involving  if  there  are    a  large 
number  of  state  and  control  variables  present  in  the  system. 
Fuzzy  interpolation  is  dealt  with  more  detail  in  Cll. 


5.2.2.  Method  of  Kacprzyk  -  Branch-And-Bound  Method 

Kacprzyk  proposed  a  much  simpler  method  for  fuzzy  system 

under  control  using  branch-and-bound  method  in  CIO].   The 

method  i s  as  follows: 

Consider  the  same  fuzzy  automata  described  in  the  previous 

section  where  the  fuzzy  transition  equation  is 

X    =  F(X  ,  U  ) 
t+i       i'   t 
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where    X   ,X  c    X    =    Cx    ,  .  .  .  ,  x    >     is    the    -fuzzy    state    at    time    t 

and    t  +  1    respectively    and    u_    e   U    =    tu   ,  .  .  .  ,  u    >     i  s    a    non--fuzzy 

tin 

control,  is  a  particular  case  o^^  the  fuzzy  system  under 
control.   The  temporal  evolution  o-f  the  fuzzy  state  at  the  end 
of  stage  t  is  given  by 

At  each  control  stage  t,  the  control  u  is  subjected  to  a 

fuzzy  constraint  Pj,<x^)  is  imposed  on  the  final  state  x  c  X. 

The  initial  fuzzy  state  i s  x   e  X,   The  branch-and-bound 

o 

technique  works  the  way  it  is  described  in  section  4.2.4. 

The  final  fuzzy  state  x   can  not  be  introduced  directlv 

T  ' 

into  Pg<x^)  as  the  constraints  are    of  point  values  and  both 
the  final  state  and  goal  are  fuzzy  sets  (remember  that  the 
control  is  non-fuzzy).   Therefore  some  trickery  is  needed  to 
apply  the  general  framework  in  case  of  this  fuzzy  system.   It 
is  essential  to  attain  the  final  fuzzy  state  x   as  close  as 

T 

possible  to  the  fuzzy  goal  g''.   A  measure  of  closeness 
(similarity)  between  x^  and  g'"  does  in  fact  serve  the  same 
purpose  as  Pjj<x^)  for  a  non-fuzzy  final  state  x  ,  provided 
that  it  takes  on  a  value  from  C0,1]  with  1  indicating  that  x 

T 

and  G  are    the  same,  zero  indicating  that  they  are  fully 
different,  and  intermediate  values  indicating  varying  degrees 
of  differences. 

Kacprzyk  proposed  normalized  linear  (Hamming)  distance 
which  is 
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d^<x^,G''>  =  1/n  2:  AJ^  (s.)  -  ;^'^(5^)  (5.12) 

T 

where  s.  are    the  states. 

x. 

For  the  normalized  distance  the  following  relation  holds 
good:   if  «   =  G^,  then  d (x  ,G^)  =  0,  and  the  more  x   differs 

T  T  T 

from  6  ,  the  closer  d(x  ,G^)  to  1.   If  ^^    (x  )  is  the  measure 

T  02    T 

of  closeness  of  final  fuzzy  state  to  fuzzy  goal,  then 

Ai^^^x^)  =  1  -  d(x^,G'')  (5.13) 

Then,  the  fuzzy  decision  becomes 
/Li(u,...u    lx)=    max     (l/**(u  )  ^  ...^   u^~*(u   )  ^ 

D    O        T-1  '   O  CO  c      T-1 

u  ,  .  .  .  ,  u 

O       '   T-1 

^  t/    (x  )  (5. 14) 

az      T 

The  case  of  fuzzy  control  may  be  dealt  with  by  the 
approximation  of  the  fuzzy  controls  applied  by  the  reference 
fuzzy  controls,  so  as  a  finite  number  of  possible  controls 
exist  which  is  a  pre-requi si te  to  construct  and  handle  the 
decision  tree. 


5.3.  Difficulties  in  Both  the  Approaches 

Both  approaches,  Baldwin  and  Pilsworth's  and  Kacprzyk's, 
suffer  from  a  difficulty.   The  basis  of  both  the  methods  is 
the  fuzzy  state  mapping.   To  obtain  this  mapping,  one  has  to 
know  all  possible  states  the  system  goes  through  C23].   This 
is  hardly  the  case  in  actual  practice.   In  most  cases,  the 
decision  maker  is  not  aware  of  all  possible  states.   In  most 
of  the  backward  flow  problems,  for  example,  resource 
allocation,  travelling  salesman,  and  network  problems,  it  is 
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possible  to  enumerate  all  possible  states.   But  it  is 
extermely  difficult  and  often  impossible  to  enumerate  all 
possible  states  in  most  of  the  forward  flow  (progressive  in 
time)  problems  with  the  exception  of  a  few  simple  cases  of 
inventory  problems.   Also,  solving  the  problem  by  both   the 
approaches  becomes  tedious  if  the  number  of  all  possible 
states  is  extremely  large.   In  such  cases,  one  has  to  look  for 
a  more  intuitive  and  simpler  approaches  in  order  to  solve  the 
problem  of  fuzzy  system  under  control  using  Bel  1 man-Zadeh 
model.   The  problem  which  we  are  currently  solving  falls  into 
this  category. 

5.4.  Solution  to  the  Problem  -  Fuzzy  System  Under  Control 

In  this  approach,  for  the  sake  of  simplicity,  we  assume 
non-fuzzy  control.   In  order  to  eliminate  unnecessary 
computations,  we  use  the  branch-and-bound  method  with  some 
modifications  as  the  original  method  described  in  the  previous 
section  can  not  be  applied  directly  to  solve  our  problem.   The 
major  features  of  our  problem  in  the  framework  of  fuzzy  system 
under  control  are: 

1.  The  percentage  production  factors  of  alternative  components 
are  fuzzy  sets.   This  results  in  fuzzy  initial  state. 

2.  All  constraints  except  that  of  machine  utilization 
(overall)  are  imposed  on  state  variables.   The  constraint  of 

overall  machine  utilization  is  applied  on  the  control. 
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3.  The  re-ference  sets  for  overall  goal  and  constraints  are 
same  as  in  the  section  4.2.3. 

4.  As  we  shall  see  later,  the  state  variables,  after  imposing 
constraints,  yield  type  2  fuzzy  sets  -  which  are  difficult  to 
handle 

5.  Since  we  have  assumed  non-fuzzy  control,  the  grades  of 
membership  of  constraints  imposed  on  control  (overall  machine 
utilization)  have  point  values  whereas  the  state  variables  are 
fuzzy  and  do  not  have  point  values  of  their  grades  of 
membership . 
Given: 

Two  machines  with  production  capacity  of  600  per  period 
each  and  a  minimum  batch  quantity  of  400  each. 

Production  plan: 

Period  Quantity 

1  1000 

2  900 

Demand:  CI  -  15*/.,  C2  -  657.,  C3  -  207.,  non -fuzzy. 
Period  Quantity 

1  1000 

2  900 
Production  alternatives  for  the  period  1 

u         u         u 

12  3 

Machine  1  400  500  600 
Machine  2  600  500  400 
Percentage  production  factors  for  alternative  components: 
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5.   = 
V  1 


Machine 

1: 

Compo 

P 

1 

nent 

CI 

15 

C2 

55 

C3 

30 

Machine 

2: 

P 

1 

Compo 

nent 

CI 

10 

C2 

70 

C3 

20 

0.7       0.3 

20  15 
60  60 
20        25 

0.7        0.3 

15         20 

s   =        C2  70       60         65 

25         15 

Note  that  for  the  sake  of  simplicity  we  have  chosen  same 

membership  -function  values  -for  all  alternative  components. 

This  would  make  the  computations  easier. 

Cal cul ati  ons: 

Stage  0,  u  =1: 
1 

Total  quantity  to  be  manufactured  on  machine  1  =  400 
Total  quantity  to  be  manufactured  on  machine  2  =  600 
Now,  the  quantity  of  CI,  C2,  and  C3  manufactured  on  machine  1 


=  Q.  X  s,   =  q. 


0.7      0.3 


comp 

=  400  X  s    =           CI 

60 

80 

60 

C2 

220 

240 

240 

C3 

120 

BO 

100 

Similarly,  for  machine  2, 

we 

h 

ave 
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CI  60       90     120 

^iii  ^    ^^^    ^   ^u  ^   "-^                 ^20      360  390 

C3                 120      150  90 

Total  quantity  of  each  alternative  components  made  =  Quantity 
made  on  machine  1  +  quantity  made  on  machine  2 


Comp 

(J 

1 

0.7 

0.3 

CI 

120 

170 

180 

C2 

640 

600 

630 

C3 

240 

230 

190 

which  is  equivalent  in  fuzzy  notation, 
q^^^  =  1/120  +  0.7/170  +  0.3/180 
q     =  1/640  +  0.7/600  +  0.3/630 

2  11 

q^    =  1/240  +  0.7/230  +  0.3/190 

3  11 

Similarly,  the  total  quantity  of  each  alternative  component 
made  for  u   and  u  are: 

2  3 

For  u  : 

1 

q^j2    ^    1/125  +    0.7/175    +    0.3/175 

q^j^    ^    l/<^25  +    0.7/600    +    0.3/625 

Qg^^  =  1/250  +  0.7/225  +  0.3/200 
For    u   : 

3 

q^^g  =  1/130  +  0.7/180  +  0.3/170 
q^^^  =  1/610  +  0.7/600  +  0.3/620 
qg^g  =  1/260  +  0.7/220  +  0.3/210 


Now,  subtracting  the  demand  (CI:  150,  C2:  650,  C3:  200>  fr 
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om 


the    total     quantities,     we    get    inventory    at    the    end    of    period     1, 

2 

i.e.  ,     I         =   I  q...    -   d.  . 

tlk  tjk  vj 

■for    u  : 

CI:    1/-30  +  0.7/25  +  0.3/30 

I         =       C2:  1/-10    +    0.7/-50    +    0. 3/-20 

111 

C3:  1/40    +    0.7/30    +    0.3/-10 

CI:    1/-25  +  0.7/25  +  0.3/25 
C2:    1/-25  +  0.7/-50  +  0.3/30 
C3:    1/50  +  0.7/25  +  0.3/0 


■for    u 


L12 


•for    u 


I 


x.1.3 


CI: 
=      C2: 
C3: 


1/-20    +    0.7/30    +    0.3/20 
1/-40    +    0.7/-50    +    0.3/30 

1/60    +    0.7/20    +    0.3/10 
Now,     using    the    reference    graphs    4.6,     4.7,     4.8    for    constraint     1 
(minimize    shortages),    we   get 

for    u   : 

1 

(J°     (u°)     =    1/0.4    +    0.7/1    +    0.3/1 
cii       1 

/J**      (u°)     =     1/1     +    0.7/0.82    +    0.3/1 

C12         1 

^J°       (u°)    =    1/1    +    0.7/0.82    +    0.3/1 

C13         1 


f  or    u   : 

2 

p        (u    ) 

Cll         2 

u         <u     ) 

C12         2 

/J  <u     ) 

C13         2 


1/0.67  -t-  0.7/1  +  0.3/1 
1/1  +  0.7/0.82  +  0.3/1 
1/1  +  0.7/1  +  0.3/1 
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■for  u  : 

3 

M°  (u°)  =  1/0.67  +  0.7/1  +  0.3/1 

Cll         3 

Ai°      <u°)     =    1/0.92    +    0.7/0.B2    +    0.3/1 

C12         3 

/J°      (u°)     =    1/1     +    0.7/1    +    0.3/1 

C13         3 

Note    that    the    above    -fuzzy    sets    are    o-f    type    2 

Similarly,     using    the    reference    graph    4.9    -for    constraint    2 

(maximize    the    quantity    o-f    C2)  ,     we    get 

Ij"     (u°)     =    1/0.6    +    0.7/0.5    +    0.3/0.58 

C2  1 

/LJ°    <u°)     =    1/0.56    +    0.7/0.5    +    0.3/0.58 

C2         2 

p°    (u°)     =    1/0.53    +    0.7/0.5    +    0.3/0.55 

C2        3 

Now,  -for  the  first  stage, 

u  (u  )  =   max   (u   ^  u   /^  u  ^  u      ) 

D  Cll     C12     C13     C2 

u  ,  u  ,  u 

1    2    3 

P°(u'*)  =  (1/0.4  -^  0.7/1  -♦-  0.3/1)  ^    (1/1  -t-  0.7/0.82  -•■  0.3/1)  ^ 

D    1 

^  (1/1  -^  0.7/0.82  +  0.3/1)  A  (1/0.6  -^  0.7/0.5  -t-  0.3/0.58). 

Here,  we  use  the  extension  principle  for  min  operator  on  type 

2  fuzzy  sets  (already  described  in  3.2.7.).   A  sample 

calculation    of    such    operation    is    given    below: 

(1/0.4    +    0.7/1    +    0.3/1)    ^    (1/1    -1-    0.7/0.82    +    0.3/1)    = 

(1    ^    1) / (0.4    ^    1)     =    1/0.4 

(1    ^   0.7)/(0.4    ^   0.82)     =    0.7/0.4 

(1    ^   0.3>/(0.4    A    1)     =    0.3/0.4 

<0,7    A    1)  /  (1    /N    1)     =    0.7/1 

(0.7    A   0.7)/(l    A   0.82)     =    0.7/0.82 

(0.7    A   0.3)/ ( 1    A    1)     =    0.3/1 

(0.3    A    1)/  (1    /N    1)     =    0.3/1 
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(0.3  A  0.7)/(l  A  0.B2)  =  0.3/0.82 

(0.3  A  0.3)/(l  A  1)  =  0.3/1 

Now,  for  1,  we  have  max  <0.4}  =  0.4  ^  1/0.4 

for    0.7    max  {0.4,  1,  0.82>  =  1  =»  0.7/1 

for    0.3,   max  {0.4,  1,  1,  0.82,  1  >  =  1  =»  0.3/1 

.-.  (1/0.4  +  0.7/1  +    0,3/1)  A  (1/1  +  0.7/0.82  +  0.3/1)  = 

(1/0.4  +  0.7/1  +  0.3/1) 
Continuing  the  same  process,  we  get 
M°<u^)  =  (1/0.4  +  0.7/0.6  +  0.3/0.6). 
Si  mi  1  ar  1  y, 

Mj^(u^)  =  (1/0.56  +  0,7/0.56  +  0.3/0.56) 
M°(u°)  =  (1/0.53  +  0.7/0.53  +  0.3/0.55) 
See  fig.  5.1. 

Now,  we  have  to  choose  the  maximum  value  of    /j'    among  three 
alternatives  in  order  to  branch  off  to  the  next  stage.   Since 
all  the  values  of  /j^  are    type  2  fuzzy  sets,  it  is  not  possible 
to  determine  which  one  has  the  greatest  value.   In  order  to 
overcome  this  difficulty,  we  propose  a  simple  and  intuitive 
method  to  evaluate  the  greatest  value  of  three  type  2  fuzzy 
sets. 

In  order  to  compare  the  type  2  fuzzy  sets,  we  need  to 
convert  them  into  a  point  value.   This  could  act  as  an  index 
based  on  which  we  can  compare  the  type  2  fuzzy  sets.   This 
index  can  be  obtained  by  taking  a  weighted  average  of  the  two 
grades  of  memberships.   Obviously,  the  weighted  average  lies 
in  the  interval  C0,1].   A  value  of  1  indicates  total  crispness 


75 


and  the  other  values  within  the  interval  indicate  varying 
degrees  o-f  cr ispness.   Therefore,  we  can  choose  that  set  as 
the  best  which  has  the  highest  value  o-f  crispness  index. 

Consider  a  -fuzzy  set  of  type  2  as  given  below. 

A  =  C(y^,Ai(y^>  )  ,  W^,  A^iy^)  >  ,  .  .  .  ,  (y^,A^(y^))> 
where  Y^,  y^ y^,  f^^y^>  ,  •  -  •  ,  f^<'y^>     €  CO,  n,  and 

yeY,  y=/(x),  x  eX. 

n 
Then,  crispness  index  =  j'  =  1/n  J]  y./J'y)  (5.15) 

<,  =  i 

Now,  using  the  above  definition,  we  get 

^•(u^)  =  1/3  (1  X  0.4  +  0.7  X  0.6  +  0.3  x  0.6)  =  0.3333 

Simil ar 1 y, 

^(u")  =  0.3733  and 

2 

y(u**)  =  0.3553 

3 

Since  ^(u^)  =  0.3733  is  the  greatest  of  the  three  values,  we 

.       o         o#    ,  o, 

choose  u  .  I.e.  u    =  (u  >. 

Stage  2. 

Production  plan  for  the  period  2  is  900.   Using  the 
equations  2.8  and  2.9  ,  we  obtain  the  following  production 
al ternati  ves. 

Machine  Alternatives 

u        u        u        u 

12  3  4 

1  300      400      500      600 

2  600      500      400      300 
Demand  for  the  period  2  is  900  units,  i.e. 
d   :  135,  d   :  585,  d   :  180 

12  22  32 

Using  the  percentage  production  factors  for  alternative 
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components,  we  get  the  following  total  quantities  of  each 
alternative  components  produced  by  two  machines. 

for  u 

t 

q     =  1/105  +  0.7/150  +  0.3/165 

12  1 

q     =  1/585  +  0.7/540  +  0.3/570 

221 

q     =  1/210  +  0.7/540  +  0.3/570 

32  1 

for  u 

z 

q     =  1/110  +  0.7/155  +  0.3/160 

122 

q     =  1/570  +  0.7/540  +  0.3/565 

222 

q     =  1/220  +  0.7/205  +  0.3/175 

322 


for  u 


q     =  1/115  +  0.7/160  +  0.3/155 

123 

q     =  1/555  +  0.7/540  +  0.3/560 

223 

q     =  1/230  +  0.7/200  +  0.3/185 

323 


for  u 


q     =  1/120  +  0.7/165  +  0.3/150 

12* 

q     =  1/540  +  0.7/540  +  0.3/555 

224 

q     =  1/240  +  0.7/195  +  0.3/195 

324 

Now,  assuming  that  there  is  no  backlogue  from  the  period  1 
(i.e.  lost  sales),  and  using  the  equation  2.   i.e. 
^  ,1,  ^  ^,1,  "*■  ^  .,,  "  cl-,  where  if  I  ,  <  0,  I.  ,   =  0, 

>.2k       tZk       vlk       v2  ».lk       '    nk 

we  get 
for  u 

*)  =  1/-30  +  0.7/40  +  0.3/55 
>  =  1/0  +  0.7/-45  +  0.3/10 
=  1/aO  +  0.7/55  +  0.3/-15 


I    (u 

12  1 

I      (U 
22  1 


I     <U 
321    1 
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I  (u*)     =    1/-25    +    0.7/45    +    0.3/50 

122         Z 

I         (u*)     =    1/-15    +    0.7/-45    +    0.3/5 

222         2 

I        (u*)    =    1/90    +    0.7/50    +    0.3/-5 

322         2 


I        (u*)     =    1/-20    +    0.7/50    +    0.3/45 

123        3 

I        <u*>     =    1/-30    +    0.7/-45    +    0.3/0 

223        3 

I    <u*)  =  1/100  +  0.7/45  +  0.3/5 

323    3 


I   <u*)  =  1/-15  +  0.7/55  +  0.3/40 

124        4 

I         <u*)     =    1/-45    +    0.7/-45    +    0.3/-5 

224         4 

I         (u*)     =    1/110    +    0.7/40    +    0.3/15 

324         4 


Using    the   reference   graphs    -for    constraints,    we   get 
/J*      (u*)     =    1/0.32    +    0.7/1    +    0.3/1 

Cll        1 

/J*      (u*)     =    1/1    +    0.7/0.B3    +    0.3/1 

C12         1 

p*      (u*)     =    1/1    +    0.7/1    +    0.3/0.78 

CIS        1 


^J^      (u*)  =    1/0.46    +    0.7/1    +    0.3/1 

p*      (u*>  =    1/1    +    0.7/0.89    +    0.3/1 

C12        2 

/J*^g<u*)  =    1/1    +    0.7/1    +    0.3/1 


l^t.    ^^i^  =  1/0.61  -t-  0.7/1  +  0.3/1 

Cll    3 

A^^^jj<u*)  =  1/0.99  +  0.7/0.83  +  0,3/1 

/J*^3<u*>  =  1/1  +  0.7/0.83  +  0.3/1 


AJ*^^(u*)  =  1/0.76  +  0.7/1  -t-  0.3/1 
/J*^2^^^'  ""  1/0.83  +  0.7/0.83  +  0.3/1 
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M*  „<u*)     =    1/1    +    0.7/1    +    0.3/1 

C13        4 

Constraint    2    (maximize    quantity    of    C2)    becomes 
P^j,^*^!'     "^    1/0.63    +    0.7/0.5    +    0.3/0.5B 
P*j,<u*)     =    1/0. 58    +    0.7/0.5    +    0.3/0.57 
Ai*    <u*)     =    1/0.54    +    0.7/0.5    +    0.3/0.56 

CZ    3 

AJ*^^'-'!^  =  1/0.5  +  0.7/0.5  +  0.3/0.54 

Overall  machine  utilization  constraint  becomes 

^ca/^^  =  0.86;  ^^^^iu^)    =    1 

^ca/^^  =  ^5  ^caz^"^  =  ^-^^^ 

^C9i<^>  =  '■'     ^C32<"9>  =  ^-^S 

Note  that  the  grades  o-f  membership  of  overall  machine 

utilization  are  point  values  since  we  have  assumed  non-fuzzy 

demand  and  non-fuzzy  control. 

Goal  : 

Z(u°)  =  1/375  +  0.7/325  +  0.3/365  for  stage  0. 

for  stage  1,  we  have 

Z(u*>  =   1/375  +  0.7/294  +  0.7/339 


Z<u*)  =  1/350  +  0.7/293  +  0.3/ 


332 


Z(u*)  =  1/325  +  0.7/292  +  0.3/325 


Z(u  )  =  1/329  +  0.7/291  +  0.3/318 

4 

Overall  Z  values  are 

Z(u^,u*)  =  1/750  +  0.7/619  +  0.3/704 

Z(u°,u*)  =  1/725  +  0.7/618  +  0.3/697 

Z(u^,u*)  =  1/700  +  0.7/617  +  0.3/690 

Z<u°,u*)  =  1/704  +  0.7/616  +  0.3/683 
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After  using  the  reference  graph  for  Z  value,  we  have 

Li*<u°,u*)  =  1/0.75  +  0.7/0.1  +  0.3/0.52 

AjNu°,u*)  =  1/0.63  +  0.7/0.09  +  0.3/0.49 
a      22 

/j*(u°,u*)  =  1/0.5  +  0.7/0.09  +  0.3/0.45 

0    2    3 

/uNu'^.u*)  =  1/0.52  +  0.7/0.08  +  0.3/0.42 

0    2    4 

Now, 

zj  <u  ,u  )  =  u  <u  )  ^         max      (u       ^   u        ^  u       ^  u      ^   u       ^ 

D  D         2  Cll  C12  C13  C2  C31 

1111 

u    ,  u    ,  u   ,  u 

1         2         3         4 

C32  a 

p*<u**,u*)     =     (1/0.32    +    0.7/1    +    0.3/1)     A    (i/i    +    0.7/0.83    +0.3/1) 

D         2         1 

A    (1/1    +    0.7/1    +    0.3/0.78)     A    (1/0.63    +    0.7/0.5    +    0.3/0.58)     A 
A  0.86    A    1    A    (1/0.75    +    0.7/0.1    +    0.3/0.52) 

Since  the  type  2  fuzzy  sets  follow  the  associative  law,  we 
can  put  the  point  values  in  the  end,  i.e. 
•^*<u°,u*)  =  (1/0.32  +  0.7/1  +  0.3/1)  A  (1/1  +  0.7/0.83  +0.3/1) 

D    2    1 

A  (1/1  +  0.7/1  +  0.3/0.78)  A  (1/0.63  +  0.7/0.5  +  0.3/0.58)  A 

A  (1/0.75  +  0.7/0.1  +  0.3/0.52)  A  0.86  A  1 

p*(u°,u*)  =  (1/0.32  +  0.7/0.63  +  0.3/0.63)  A  0,86  A  1 

D    2    1 

Si  mi  1  ar  1  y  , 

p*(u°,u*)  =  (1/0.46  +  0.7/0.58  +  0.3/0.58)  A  1  A  0.89 

D    2    2 

p*(u°,u*)  =  (1/0.5  +  0.7/0.5  +  0.3/0.5)  A  1  Ao.75 

D    2    3 

p*(u°,u*)  =  (1/0.5  +  0.7/0.5  +  0.3/0.52)  A  1  A  0.6I 

D    2    4 

Now  we  calculate  the  crispness  index  for  the  type  2  fuzzy  sets 

to  get  a  point  value,  i.e. 

AJ*(u*',u*)  =  0.3167  A  0.86  A  1  =  0.3167 

D    2    1 

p*(u°,u*)  =  0.3467  A  1  A  0.89  =  0.3467 

D    2    2 
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/Lj*(u°,u*>  =  0.3333  ^  1  ^   0.75  =  0.3333 

D    2    3 

AJ*<u°,u*)  =  0.3353  ^  1  ^   0.61  =  0.3353 

D    2    •• 

1    O    1 

Since  /J  (u  ,u  )  has  the  highest  value  with  0.3467,  we  choose 

D    2    2 

u      for    period    2.       There-fore    the    -final     solution    is 

Stage    1 
u 

2 

Machine  1:  400 
Machine  2:  500 


Stage  0 

u 

2 

Machine 

1:  500 

Machine 

2:  500 

See  Fig. 

5.2. 

ei 


v^=Cl/.  4  +  .  7/.  6+.  4/.  6D   CI/.  56  +  .  7/.  56  +  .  3/.  563 


CI/.  53  +  .  7/.  S3  +  .  3/.  55:) 


Fig.  5.1.  Decision  tree  f 


or  stage  O. 


V  =0. 3333 
o 


V  =  0.  3167 


O. 3457        0. 3333 
Fig.  5.2  Decision  tree  for  the  exampli 


0. 3533 


O. 3353 
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CHAPTER  6 
SUMMARY  AND  CONCLUSIONS 

6. 1 .  Summary 

Alternative  components  are  very  common  in  manu-f  actur  i  ng 
environment.   They  components  arise  out  of  lack  of  total 
control  over  manufacturing  processes  and  have  same  basic 
dimensions  and  tolerance  but  different  upper  and  lower  limits. 
Since  alternative  components  form  many  options  for  a  single 
part,  it  is  difficult  to  plan  and  control  the  production  of 
such  components.   Therefore  it  is  desirable  to  maximize  the 
production  of  the  main  component  which  is  most  desirable, 
whenever  alternatives  exist,  by  choosing  a  proper  one.   In 
situations  where  a  part  is  manufactured  on  many  similar 
machines  each  having  its  own  degree  of  control,  it  is  possible 
to  have  different  alternative  plans  by  splitting  the  order 
quantity  among  the  machines  in  such  a  way  as  to  maximize  the 
quantity  of  the  main  component  produced.   Since  this  decision 
is  an  ongoing  process  over  a  length  of  a  finite  number  of 
periods,  it  becomes  a  multi-stage  decision-making  problem 
where  dynamic  programming  is  extensively  used. 

This  study  develops  a  mathematical  model  of  the  problem 
of  alternative  components  using  the  dynamic  programming 
approach.   The  problem  is  treated  under  three  categories 
namel y: 
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-deterministic  system  under  control  in  a  crisp 

envi  ronment . 

-deterministic  system  under  control  in  a  -fuzzy 

environment. 

—fuzzy  system  under  control  in  a  fuzzy  environment. 
Mathematical  models  are  -formulated  for  each  category  and  an 
example  of  the  problem  is  solved  in  each  case.   Bellman  and 
Zadeh's  model  of  decision-making  in  a  fuzzy  environment  is 
used  as  a  basic  framework  to  study  the  problem.   In  case  of 
fuzzy  system  under  control  in  fuzzy  environment,  the  existing 
approaches  ars    evaluated  in  view  of  the  applicability  to  our 
problem.   It  was  found  that  the  existing  two  methods,  Baldwin 
and  Pilsworth's  method  and  branch-and-bound  method  of 
Kacprzyk,  can  not  be  applied  to  our  problem  in  fullness..   The 
latter  method  is  used  to  solve  our  problem  with  some 
modi  f  icati  ons. 
6.2.  Discussions 

In  deterministic  system,  it  is  very  essential  that  all 
the  parameters  be  crisp.   As  we  have  seen  in  chapter  2, 
conflicting  goals  and  constraints  result  in  a  deadlock 
situation  where  the  decision  maker  is  left  with  no 
alternatives  to  choose  from.   In  most  of  the  practical 
situations,  physical  quantities,  goals  and  constraints  can  not 
be  specified  detrmi ni sti cal 1 y .   In  many  situations  there  is 
vagueness  associated  with  many  parameters  and  variables  and 
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the  probabilistic  models  can  not  handle  fuzziness  as  they  are 
designed  to  take  care  of  randomness. 

Defining  the  problem  in  a  fuzzy  environment  alleviates 
some  of  the  difficulties  by  allowing  the  goals  and  constraints 
to  be  defined  with  imprecision.   Imprecise  terms  like  "as  low 
as  possible"  and  "as  high  as  possible"  can  now  be  accomodated 
in  the  decision  model  by  fuzzifying  such  parameters. 
Obviously,  the  f uzzi f icati on  is  subjective  and  done  by  a 
person  who  i s  an  expert  in  that  area  or  the  one  who  has  a  good 
knowledge  of  the  syatem.   Fuzziness  of  goals  and  constraints 
can  be  represented  by  the  reference  graphs  which  have  an  upper 
limit  or  an  aspiration  level  and  a  lower  limit  or  the  worst 
case  which  are    linked  by  a  linear  relation.   The  decision 
yields  a  point  value  of  the  grade  of  membership  and  the 
alternative  which  has  the  highest  grade  among  all  the 
alternatives  is  chosen  as  a  decision. 

Since  Bellman  and  Zadeh's  model  of  multi-stage 
decision-making  in  a  fuzzy  environment  is  designed  for 
backward  flow  problems,  it  can  not  handle  forward  flow 
problems  with  high  dimensions.   Kacprzyk's  branch-and-bound 
method  offers  a  simple  approach  for  forward  flow  problems  with 
high  dimensionality.   Instead  of  doing  a  blind  traversal  of 
all  possible  alternatives,  it  chooses  the  most  promising  paths, 
thus  eliminating  the  unpromising  alternatives.   Since  our 
problem  can  have  a  high  diment i onal i ty ,  if  there  are  more  than 
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two  machines  and  the  lot  sizes  ars    smaller  instead  o-f 
multiples  o-f  100,  branch-and-bound  method  appears  to  be  the 
most  e-fficient  method. 

We  see  in  case  of  fuzzy  system  under  control  that,  as 
the  number  of  fuzzy  parameters  increase,  the  problem  gets  more 
and  more  complicated.   Both  the  methods  presented  require  a 
fuzzy  state  mapping  to  start  with.   This  means  that  the 
decision  maker  needs  to  enumerate  all  possible  states  to 
obtain  a  transition  matrix  for  both  state  and  control.   The 
main  problem  here  is  that,  in  most  practical  situations  it  is 
almost  impossible  to  enumerate  all  possible  states.   For 
example,  if  the  inventory  level  is  considered  as  a  state 
variable,  then  the  number  of  all  possible  levels  attained  by 
inventory  will  be  infinitely  large.   In  a  very  highly 
uncertain  situation,  the  inventory  level  may  range  from  10 
units  to  100,000  units.   In  such  a  case,  there  would  be  99,990 
possible  states.   Suppose  if  the  control  variables  are    10, 
then  the  mapping  matrix  would  be  of  999,900  x  999,900  1   Since 
this  is  the  case  in  most  of  the  forward  flow  problems,  it 
would  be  impossible  to  apply  both  the  methods  to  the  problem. 
In  order  to  solve  such  problems,  it  is  essential  to  look  at 
the  basic  underlying  philosophy  governing  the  method  rather 
than  applying  the  method  in  its  rigid  framew("l<.   One  has  to 
look  for  a  simple  and  often  heuristic  approach  to  obtain  a 
rational  solution  for  a  problem.   The  crispness  index 
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developed  in  chapter  5  is  one  of  such  attempts.   We  have  used 
the  basic  concept  of  the  decision  being  the  confluence  of 
goals  and  constraints  in  our  approach.   Since  the 
branch-and-bound  method  cuts  down  the  number  of  traversals 
drastically,  this  method  has  been  adopted.   We  can  see  that 
even  with  the  assumption  of  non-fuzzy  control  we  ended  up  with 
type  2  fuzzy  sets  for  which  methods  of  comparison  are  not 
available.   One  can  imagine  the  complexity  that  would  have 
arisen  had  we  used  fuzzy  control  and  fuzzy  demand.   We  can 
say  say  that  our  method  of  obtaining  the  crispness  index  by 
weighting  the  membership  grades  of  the  type  2  fuzzy  set  is 
intuitive,  appropriate  and  serves  the  purpose. 

Until  now,  only  two  methods  are    available  for  fuzzy 
system  under  control.   Both  are  applicable  for  certain 
specific  types  of  problems  where  it  is  possible  enumerate  all 
possible  states  to  form  a  state  transitional  matrix.   However, 
this  study  seems  to  be  the  first  attempt  to  solve  the  problem 
of  alternative  components.   And  the  limitations  of  the 
existing  approaches  are    evidently  clear  in  this  study.   Much 
work  needs  to  be  done  in  the  direction  of  obtaining  simpler 
and  intuitive  methods  which  can  be  applied  to  a  wide  variety 
of  problems. 
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ABSTRACT 

Alternative  components  are  very  common  in  manu-f  actur  ing 
environment.   They  are  the  result  o-f  lack  of  control  over 
manu-f  actur  i  ng  processes.   They  have  the  same  basic  dimension 
and  tolerance  but  di-f-ferent  upper  and  lower  limits.   Since 
alternative  components  form  many  options  for  a  single  part,  it 
is  difficult  to  plan  and  control  the  production  of  such 
components  and  their  mating  parts.   In  most  of  the  cases, 
achieving  complete  control  over  processes  is  impossible  or 
costs  are  prohibitive   and  the  components  Are    made  on  many 
similar  machines  each  having  its  own  degree  of  control.   Then, 
it  is  desirable  to  maximize  the  production  of  the  main 
component  by  suitably  allocating  the  order  quantities  to  the 
machines.   Since  this  decision  is  an  ongoing  process  over  a 
length  of  finite  number  of  periods,  it  becomes  a  multi-stage 
decision-making  problem  where  dynamic  programming  is 
extensively  used. 

This  study  develops  a  mathematical  model  for  the  problem 
of  alternative  components  using  the  dynamic  programming 
approach.   The  problem  is  treated  under  three  categories 
namel y: 

-deterministic  system  under  control  in  a  crisp 


envi  ronment . 

-deterministic  system  under  control  in  a  -fuzzy 

envi  ronment . 

-■fuzzy  system  under  control  in  a  fuzzy  environment. 

Mathematical  models  are  formulated  for  each  category  and 
an  example  o-f  the  problem  is  solved  in  each  case.   Bellman  and 
Zadeh  model  oi    decision-making  in  a  fuzzy  environment  is  used 
as  a  basic  -framework  to  study  the  problem.   In  case  of  fuzzy 
system  under  control  in  a  fuzzy  environment,  the  existing 
approaches  of  Baldwin  and  Pilsworth  and  Kacprzyk  are  evaluated 
in  view  of  applicability  to  our  problem.   They  were  found  not 
to  be  completely  applicable  to  our  problem  in  their  fullness 
as  the  problem  falls  outside  the  domain  assumed  by  both  the 
methods.   The  branch-and-bound  method  of  Kacprzyk  is  used  in 
solving  our  problem  with  a  modification  in  the  form  of 
crispness  index  to  evaluate  the  type  2  fuzzy  sets  resulting 
from  the  fuzzy  state. 


